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(54) Methods of characterising ventricular operation and applications thereof 



(57) New methods of characterising ventricular op.^ 
erations by measuring propagation characterises of the 
repolarisation wavef ront (the T wave) are disclosed the 
methods use new descriptions of T wave Morphology. 
Dispersion (TMD), Total Cosirt R_ to _T (TCRT) and T 
wave energy residium to quantify the wavef ront charac- 
teristics, these descriptors measure the spatial variabil- 



ity of the T wave Morphology, the vector deviations be- 
tween the.depolarisation and repolarisation wavef ronts 
and the energy, of the non-dipolar. components of the 
ECG vector respectively. TCRT also provides a respon- 
sive descriptor. for. measuring autonomic, tone As such, 
has applications for : . improved pacing and autonomic 
nervous system monitors. 
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Description 

[0001] The present inventions relate to methods of characterising ventricular operation In particular, but not exclu- 
sively, they re!ate to a system for quantifying abnormalities of an electrocardiogram and to a method and an apparatus 
5 for measuring such abnormalities, The present inventions also extend to an operating system for a computer, to a 
computer program and to media having stored thereon a computer program for putting the inventions into effect. Other 
applications include. use of the algorithms in pacemakers and heart monitors. The, inventions share a common link of 
characterising differences in the wavefront of the repolarisation wave. 

[0002]. Electrocardiographic patterns of the heart's movements have-been well. studied. Ah electrocardiogram (ECG) 
10 records the changes in electrical potential associated with the spread of depolarization and repolarisation through the 
heart muscle. In a normal healthy patient, depolarisatibn starts in an area of the right atrium called the sinoatrial node 
and spreads through the atrioventricular node and into the ventricular muscle via specialised conduction tissue. causing 
the two atria and the two ventricles to contract" During repolarisation, the atria and ventricles relax and. refill with blood 
The depolarisation of the atria is responsible for the P wave of an ECG and depblarisatipn of the ventricles results in 
is the QRS complex. Repolarisation of the atria coincides with the QRS complex so it is not seen Repolarisation of the 
ventricles, however, is. seen as the T wave 

[0003] ECG's are typically recorded using a standard arrangement of 12 leads, 6 (the I, II, III, VR, VL ; VF leads) 
looking at the heart in different directions in an approximately vertical plane of a body in an upright position and 6 (the 
VI, V2, V3, V4, V5 and V6 leads) looking at the heart in different directions in an approximately horizontal plane Using 
20 such an arrangement of leads, the spread of the waves of electrical potential associated with depolarisation and re- 
polarisation through the three dimensional space of the body can be recorded. 

[0004] The spread of these waves through the heart is often described by vectors. For example, the average direction 
of spread cf the depolarisation; wave through the ventricles as seen from the front of the body is. called the cardiac axis 
and the direction of this axis has long. been used, to indicate different abnormalities of the heart. 

25 [0005] To study abnormalities associated, with ventricular repolarisation, a number of data processing techniques 
have been proposed to measure, for exarhple, the QTintervai, i e the interval between the beginning of depolarisation 
and the end of repolarisation of the ventricles Interlead variability of the QT interval durations in standard 12 lead ECG 
recordings has also been studied. However, whilst these measurements may provide some diagnostic assistance, 
concerns have been raised abbut the poor reproducibility of results 

30 [0006] Studies have also tried to quantify the inhomogeneities in the ventricular repolarisation patterns by evaluating 
the complexity of the T wave morphology using eigenvalues associated with the principal components of ECG, meas- 
ured over a period of 24 hours The direction of the ECG vector during T wave in the 3D.physicai (x,y,z) has also been 
shown to have some predictive value 

[0007] However, there is still a need for further measurements which may provide a more accurate technique for 
35 identifying certain conditions, particularly those which affect repolarisation of the ventricles. A problem with known 
methods, for example, is that they only quantify global variations in the T wave ratherthanspatial variations in individual 
waves, that is the synchronicity of the T wave, as observed from different locations on the body, is not observed 
[0008] Thus, viewed from a first broad aspect, a first invention described herein provides a method of quantifying 
abnormalities of an electrocardiogram observing repolarisation patterns from different locations on a body, wherein the 
46 abnormalities are quantified by a measure o1 the synchronicity of the repolarisation patterns as. observed from those 
different.focationson the body In other words, this is a measure of the homogenity. of the spreadof repolarisation waves. 
[0009] Unlike depolarisation of the heart muscle, the repolarisation of individual ceils is not triggered by neighbouring 
cells but is instead a time dependent process. If repolarisation patterns, asobserved in different locations on the body, 
lack synchronicity, then this can be indicative of certain heart complications 
45 [001 0] By quantifying these abnormalities, it may be possible to use the data to assist with diagnosis or to identify 
patients at most risk or classify them into different categories of risk. This may be pi great importance in determining 
whether certain treatments should be offered to a patient, for example The data, could also be used to trigger an alert 
in a monitoring device 

[0011] The homogenity of the spread of repolarisation vyaves can be measured by quantifying the spatial variability 
so of the ventncular repolarisation patterns i e the spatial variability of the T wave 

[0012] Thus viewed from a second aspect, the first invention provides a method of characterising ventricular oper- 
ation, comprising the steps of: 

recording a signal monitoring the propagation of a repolarisation wave, 
55 determining a vector which is representative o1 the wavefront of the. repolarisation wave, and 

determining a measure of the spatial variation of the repolarisation wavefront. 

[0013] In one preferred embodiment, it provideed a method of quantifying abnormalities of ventricular repolarisation 
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by determining a measure of the spatial T wave morphology variation; ; ' 

[0014] -Preferably the spatial T Wave morphology variation is quantified by measuring the T wave -Morphology .Dis- 
persion 

[0015] Preferably this is achieved by determining vectors describing the contributions which the signals from each 
5 ■ lead (often referred to as the channels of an EGG) makes to the T wave. The angles between these vectors are then 
calculated arid a mean value is deterniihed: This mean vklu'e of the/angles provides a measure of the. spatial T wave 
morphology variation.- The smallerthe value, the closer the T wave morphologies will.be in the signals.of. the individual 
• leads.. V ' : - = • •• : '/ "V;. : ' •• ' •*" ■ ".. ' ■ • v'V." : ; . '■ 

• [0016]. Preferably , the ECG signal- is morphologically filtered to improve the/signal- to noise .ratio. . In one preferred 

io. . embc^irrient/this consists of the steps of decomposing the T wave using a technique such as Singular Value Decom- 
position, filtering by keeping only, the two most significant signal components, and applying a DC compensation A 
. preferred DC compensation is provided by subtracting an average of the start and end signal components during the 
QRS complex and T wave The morphologically filtered T wave is then preferably rescaled to equalise energies in the 
different component directions. The corresponding reconstruction parameters are calculated to determine the vector 

is contributions of each of the. ECG leads. The angles between each pair of the vector contributions is then calculated 
and the mean determined Most preferably the contribution in. respect of lead VI is ignored because the T wave mor- 
phology in this lead is generally different than that of other channels, irrespective of any. clinical background, mainly 

. .. due to the position of the VI. electrode, and by ignoringjhis. component, it has ^e effect of enhancing the predictive 
value of the T wave morphology dispersion descriptor 

;2Q [0017] The main reason for initially decomposing the data matrix is to find an optimum representation of the ECG 
signals upon which the measurements, can be performed In this way, the system does not use the standard XYZ axes 

• ■* of the.body but finds an optimally-constructed orthogonal system to represent the 1 2 lead ECG In a preferred embod- 

: ». iment, therefore, the first invention can be seen a$ providing a method for looking £t the vector. representation of each 
of the standard electrocardiographic leads in an optimum dimensional, vector space in which the ECG signals can be 

25 represented ; and comparing the angles between the vectors of individual standard leads. 

[0018] The spatial T wave. morphology variation may provide.a useful descriptor when .it is determined for the whole 
of the T wave, the first half of the T waye r the second half .of the T wave or any other portion or combination of portions 
of the T wave 

[0019] The present inventions are not limited to standard 12 lead .electrocardiograms, although this is preferred, but 
30 extend to electrocardiograms produced from only three or more leads In certain applications, it may be useful to use 
the electrodes of a pacemaker to record an electrocardiogram signal. In such situations, the positions and numbers of 
the electrodes would hot usually correspond with the. arrangement pf standard leads, for example Whilst, it is preferred 
to view the waves in three dimensions, because, research up to date suggests that approximately 99% Of 12-lead ECG 
energy can be represented in a 3D space, the inventions are applicable to situations where the heart is. viewed in any 
35 dirhehsions, greater than or equal to two ; t _ 

[0020] In a conventional 12 lead ECG, only 8 (!, .11, WV2, V3, V4..V5 and V6) of the signals are independent The 
other 4 signals (III, VR, VL, VF} are algebraically dependent on the other leads so, rf desired, may. be generated by 
data processing methods rather than measured as such As explained above, it is most -preferred to use signals only 
from leads I, II, V2, V3, V4, V5 and V6, and to ignore the signal ffomVI in order to concentrate the abnormalities seen 
40 m the T wave ' The position of the leads, although having an effect on the value produced by the. descriptor?, is not 
critical to the inventions; Whilst the inventions may. be described with reference to the standard ECG leads, this- is not 
intended to limit the inventions to just those positions These positions are preferred, however,- since they have become 
"standard" measuring points throughout the world. 

..[0021] It has also been found that, comparing the spread of depolarisation through the ventricles with the spread of 
45 repolarisation can provide useful information 

[0022] Thus, in accordance with a first aspect of a second invention described herein, there is provided a method of 
characterising ventricular operation, coriprising the steps of: 

recording a signal monitoring the propagation of depolarisation and. repolarisation waves, 
so determining vectors which are representative of the direction. of the wavefronts of the depolarisation and repolari- 

sation waves, and 

determining a measure of the deviation between those vectors. 

[0023] Thus, the present invention can be seen to provide a method of quantifying abnormalities of an electrocardi- 
55 ogram observing the spread of depolarisation and repolarisation waves through the ventricles, wherein the abnormal- 
ities are quantified by comparing a property of the depolarisation wave with a property of the repolarisation wave where 
preferably the abnormalities are quantified by a measure of the vector deviation between the ventricular depolarisation 
and the ventricular repolarisation waves. 
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[0024] Described in other terms, a method of the second invention may compare the direction of the depblarisation 
wave (i e , the QRS part of the ECG) with the repolarisation wave (i.e the T wave) This may be achieved by cqmparthg 
the angles between principal vectors of the ventricular depolarisation and repolarisation waves ; comparing.the angles 
between a principal vector of the ventricular depolarisation wave and the ECG. vectors during ventricular repolarisation, 
comparing the angles between the ECG vectors during ventricular depolarisation and a principal vector of the ventricular 
repolarisation wave, or comparing the angles between the ECG vectors during ventricular depolarisation with those 
during ventricular repolarisation. The angles may be compared for the whole of a wave or just a portion of a wave.or 
any combination of portions of the waves: For example, in one embodiment the angles between the depolarisation and 
repolarisation. vectors are compare^'for portions of the waves which .spanthe peak energy values, but it may be pre- 
ferred in some instances. to look at and compare other portions of the waves which would correspond to depolarisation 
and repolarisation occurring in different regions of the heart muscle. 

[0025] Preferably the vector deviations are determined using the optimally constructed representation of the ECG 
signals discussed above 

[0026] In healthy patients, the principal vectors would, generally speaking, only deviate by up to about 30° , In patients 

is having hypertropic cardiomyopathy.(HCM), for example, vector deviations greater than 90° may be seen. These de- 
viations can be distinguished over inversion of the T wave, for example, which would result in angles closer, to 180° 
[0027] Again it is preferred tp conduct certain data processing steps before the angles of the vectors describing the 
QRS complex and T wave are compared Firstly the data matrix describing the signal is decomposed, again preferably 
by Singular Value Decomposition The decomposed signal components are ranked in order of their significance in 

20 terms of the energy of the ECG vector that they represent Thus the first signal .component contains the most energy 
in a first direction. The second signal component contains the next most energy, in a second direction which is perpen- 
dicular to the first. The third signal component containsthe ;next most energy in a third direction which is. perpendicular 
to the first and second directions. Where eight of the independent ECG channels are recorded,- the .ECG vector can 
be decomposed into an eight dimensional orthogonal space. When rneasuring.the vector deviations.of depolarisation 

25 and repolarisation vectors, a good approximation may be made by only measuring the first two or three of the decom- 
posed signal components since these can account for. more than 99% of the total, energy of the 1 ! 2-lead ECG signal 
The QRS complex and T wave are localised by making use of the Variation of the instantaneous ECG energy. The 
method does not depend on accurate localisation of the QRS complex and the T wave. This method 6f detection 
provides an example of many possible ways. 

30 [0028] Both the vector representations of the QRS complex and T wave follow an approximate loop, in the constructed 
space Vectors can he determined which represent the maximum energy of the T wave and QRS complex, and the 
angles between therri compared More preferably a vector describing the maximum energy of the T wave. is compared 
to the vectors describing the QRS complex for a set of points at the peak of ventricular repolarisation (which, for 
example, can . be determinedly tracing the instantaneous ECG. energy). In the most preferred embodiment, the vector 

35 deviation between the depolarisation and repolarisation waves is measured as the average of th e cosines of the angles, 
between the vector describing the maximum energy of the T wave and the vector describing the. wavef ront of the QRS 
complex, the angles being determined in constructed space The measurement of the vector deviation in terms of the 
cosine of the angle is .referred to herein as TCRT - total cpsine R_to_T; 

[0029] The second invention introduces the idea of considering depolarisation and repolarisation processes of the 
40 heart muscle simultaneously and comparing them While it may be described as a comparison of ECG vectors observed 
during these processes, in an appropriately constructed vector space, this is not intended, to limit the invention to just 
ECG vector comparison and to processing of standard ECG leads. 

[0030] Comparing the propagation directions, of depolarisation and repolarisation has revealed some interesting de- 
tection properties In particular, TCRT has proved to be more sensitive to autonomic changes of ventricular repolari- 

45 sation than other known descnptors such as ventricular gradient and QT dispersion From investigations, it has been 
found that TCRT responds quickly to changes in the position and activity of the patient with distinct ranges or levels 
of descriptor values being obtainable for different autonomic tones This descriptor may be used. to. check for abnor- 
malrties in ventricular depolarisation and repolarisation under different autonomic conditions, thereby.providing a fuller 
picture tp assist with diagnosing defects. TCRT has been found to be useful as a predictor for mortality in patients 

so which have suffered acute myocardial infarction and as a predictor for arrythmias. The descriptor could be used in a 
monitor carried by the patient or in equipment in an intensive care unit to warn the patient or medic by means of an 
alarm when the TCRT is pushed to a danger level as a result of changes in the autonomic tone, for example, caused 
through exercise or trauma. TCRT could also be used to check that the patient has a propisrly f unctioning autonomic 
system, for example, prior to the administration of anaesthetic agents before surgery. 

55 [0031] TCRT is in effect able to provide a measure of the autonomic tone of a patient It could be used to control 
pacing of a. pacemaker making it more responsive to the patient's needs by responding to changes in the autonomic 
tone TCRT is increased by physical or emotional stress during fixed rate pacing and is decreased by an increase in. 
pacing rate TCRT could be implanted in a closed loop rate-adaptivef eedback system. At times of physical or emotional 



4 



stress TCRT would increase, triggering an increased pacing rate to decrease TCRT to resting level. - ■ 
[0032] TCRT could be used to monitor the effect of certain drugs and the way in which they effect the autonomic 
system When testing drugs which prolong the QT interval, TCRT could be used to monitor the patient and raise an 
alarm arrythmias are predicted or detected. It could be used to monitor changes in electrolyte of the body and other 
5 conduction phenonoma. It could even be used to control a drug delivery mechanism, administering certain drugs as 
heart-function abnormalities are detected or in respond 

[0033] Conditions such as ischemia, as well as most illnesses, will have, an effect on the autonomic. tone of the 
patient. TCRT could be used to assist in the prediction of ischemia or jn the monitoring of the progressionof a diseasei, 
.; for. example, in heart failure, patients, by providing.an: : .indication ; pf.the autonomic tone, as, well^ 
to' of the repolarisatioh wavefront; TCRT may be abje to observe autonomic changes .;caused by the onset of ischemia 
before St segment-changes are observable on an ECG or pain; is felt by the. patient It may also be useful in the 
monitoring of patients suffering from epilepsy, providing an early warning of heart function abnormalities The autonomic 
tone could be observed "to detect hypoxy conditions in a patient prone to fitting. 

[0034] . It should be noted that the possibilities mentioned above with reference to TCRT are not mtendedto be inferred 
?5 as limiting, the .present invention to the preferred situation wh ere TCRT is the average of the cosines of the angles 

between the vector describing the maximum energy of the T wave and the vector describing the wavefront. of the QRS 

complex For example, TCRT may compare the angles between the sets of vectors describing the depolarisation and 
. 4 ;repplarisation..wavefronts with respect to time or may compare the vectors describing the maximum. energy of depo- 

larisation and repblarisation to each other While the use of cosine, provides an effective way of separating the angles 
20 -between the vectors associated with abnormalities from those observed in normal patients, other operators may be 

used to separate the data 

[0035] As mentioned above, of the standard 1 2 lead ECG signals, .8 cire independent Thus, it is possible to. describe 
the T wave as an 8-by-n matrix M, with, each .row corresponding to a standard ECG channel (I, Jl, yi , V2, V3, V4, V5, 
V6) and n being the number of samples Performing Singular Value Decomposition ..on the matrix M generates a signal 

25 vector representing the progress of a T wave in 8 dimensions, where each dimension can be regarded as a component 
of the signal, associated with a fraction of the total energy of 1 2-ieadECG. For most purposes, as mentioned above, 
only the first two or three components are normally used, since these may account for over 99% of the total T wave 
energy in 12-lead ECG signals However it has been found that comparing the energy of the most Signjficant compo- 
nents with, the energy of the other components provides a further useful descriptor that can be used during analysis 

30 of the ECG 

[0036] Thus according to a third invention described herein, there .is provided a method of quantifying abnormalities 
of an electrocardiogram having a plurality of independent signals, in which the signals are decomposed to obtain a 
signal vector having two or more signal vector components, wherein the energy of the components is Compared. 
[0037] : Preferably, the components are arranged substantially in order of decreasing signal energy and the energy 

35 of the most significant components representing the majority of the signal energy is. compared to the energy of the 
other components. Preferably the electrocardiogram records. 8 independent signals and the signal vector has 8 com- 
ponents and the first 3 components of the signal vector, thus constructed, represents the majority, of the signal energy. 
[0038] The third invention introduces the idea to transform the ECG signals into an optimally constructed space which 
represents the ECG.snergy in a minimum, dimensional space of orthogonal components and to assess the residual 

40 energy that is left outside the three dimensional space of the first.three signal components. In a possible implementation, 
as described above, this corresponds to assessing the relative values of the 3 highest singular values of the matrix M 
to the other singular values. The first three signal components represent the dipolar components of the ECG vector 
and the remaining signal components correspond to the non-dipolar components. By measuring the power in the 
orthogonal components outside the three dimensional space, it is possible to measure the power of the non-dipolar 

& components which provides a rneasure of the local repolarisation abnormalities , However with a different set of input 
signals and/or different parts of the ECG signal, the separation/nufnber of significant components may vary: 
[0039] The inventions described above also extend to ah apparatus which, is. programmed with an algorithm. to proc- 
ess data from an ECG in accordance, with any of the.described inventive methods. The apparatus could be a computer, 
for example, programmed in a particular, way or could be a plug-in box for an ECG. apparatus or an ECG. apparatus 

50 provided with means to calculate these descriptors, and display the result to an operator. Furthermore, the inventions 
extend to an operating system or a computer program having an algorithm to process data from an ECG in accordance 
with the described methods and to media having such a computer. prog ram or operating systerh stored thereon Thus, 
Ihe inventions extend to a computer program product which is direotfy loadable in.to.the internal memory of a digital 
computer, comprising software code portions for performing the steps of the aforendescribed methods when the product 

55 is run on a computer 

[0040] The present inventions will now be described by way of example only with reference to a preferred embodiment 
and the accompanying drawings, in which. 
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| Figure 1 shows a flow chart of preferred algorithrhs; 

I Figures 2a and 2b shoWs.B input and 8 decomposed ECG signals (which can be obtained by Singular. Value 

| Decomposition), respectively; 

j Figure 3 illustrates approximate QRS and T wave.detection using ECG energy, E 3D , which is calculated from St . 

5 s 2 and S3, the most significant 3 decomposed ECG signals; 

Figure 4 provides an. example of the t loop Path of the tip of (the most significant point of the ECG vector in 
the decomposition space) on theu^ plane; . 
. Figured illustrates an.example of the cell weights D's,. assigned 1 to 100 ; egual cells in v,u 2 plane- (excluding the 
, : . zero weights), in increasing order^ showing the 30* celi.above the threshqld, the SO^ cell {the one shown in.Figure 
10 4) is .selected as the T wave end, 

Figure 6 shows the 3 most significant decomposed, time-orthoggnal, channels of a T wave, 

Figures 7a and 7b illustrate a T wave loop and reconstruction vectors of each standard ECG lead for a norma! and 

HCM patient respectively; 

Figures 8a and 8b illustrate QRS and T wave loops for a normal and HCM patient respectively, 
rs Figure 9 illustrates the reproducibilities of all the descriptors,, as measured by the ratio of individual variation to 

total variation in 1 0 supine recordings, from 76 normal and 63 HCM subjects (Grey: Normal subjects, Black. HCMs), 
Figures 10a and 10b show the variations in descriptor value for ventricular gradient angle, ventricular, gradient 
magnitude/TCRT, RR interval, maximum QT and QT dispersion for.the postural changes of resting supine position, 
followed by sitting, unsupported standing, supine and standing position, 
20 Figures 11a and 11b show the variations in descriptor value for ventricular gradient angle, ventricular gradient 

magnitude, TCRT, RR Interval, maximum QTand QT dispersion during Valsalva manoeuvre, each manoeuvre 
being, performed 3 times in supine position and 3 times in unsupported standing position, 
Figure 12 shows an implantable device positioned within the human body;. 
Figure 1 3. shows a perspective view of a preferred pacemaker, 
2s Figure 14 shows a circuit diagram suitable for operating the pacemaker of Figure 1 3, 

Figure 15 shows a perspective view of a. preferred pacemaker-cardibverter-defibnliator, 
Figure 16 shows a circuit diagram suitable for operating.the pacerhakerof Figure 16, 

Figures 17a and 1 7b show scatter diagrams of QT dispersion, values, where QTd Method 1 = range of measurable 
QT intervals, QTd Method 2 = standard deviation of measurable QT intervals, and QTd Method 3 = interquartile 
30 difference of measurable QT intervals, 

Figures 18a and 18b show QTd Method 1 and QTd Method 2 results displayed with respect to the clinical group 
of the subjects, 

Figure 19 shows the relative T wave residuum plotted with respect to the clinical group of the subject, and 
Figures 20. and 20b show the lack of relationship between QT dispersion and relative T wave residuum data 

35 

[0041] Three new approaches for the analysis of ventricular repolarisation in 12 lead electrocardiograms (ECG) will 
now be discussed in relation to a first case study: The spatial and the temporal variations of T wave morphology and 
the wavef ront direction difference between the ventricular depolarisation and repolarisation waves A minimum dimen- 
sional space, constructed by the Singular Value Decomposition of ECG. signals, is used. The spatial variation charac- 
40 terises the morphology differences between standard leads. The temporal variation measures the change of interlead 
relations during ventricular repolarisation The depolarisation and repolarisation patterns are compared using the prin- 
cipal interlead relations (wavef ront directions) that characterises them All of the descriptors are measured using the 
ECG vector in the constructed space and the singular vectors that define this space None of the descriptors requires 
time domain measurements (e g. the precise detection of the T wave offset), avoiding the inaccuracies associated with 
45 the conventional QT interval relatedparameters 

[0042] /The new descriptors have been compared with the conventional measurements provided by a commercial 
system for an automatic evaluation of QT interval and QT dispersion in digitally recorded 12 lead ECGs (QT Guard, 
Marquette Medical Systems) • The basic comparison used a set of 1100 normal ECGs The short-term intrasubject 
reproducibility of the new descriptors was compared to that of the conventional measurements in a set of 760 ECGs 
so recorded in 76 normal subjects and a set of 630 ECGs recorded in 63 patients with hyperthropic cardiomyopathy (1.0 
serial recordings in each subject of both these . sets) The discriminative power of the new and conventional parameters 
to distinguish normal and abnormal repolarisation patterns was compared using the same set 
[0043] The results showedthatthe new parameters do notcorrelate with theconventional ones (i e. assess different 
ECG qualities), are generally more reproducible than the conventional parameters based on the QT interval measure- 
rs ment, and lead to a more significant separation between normal and abnormal ECGs both univanately and in multi- 
variate regression models 

[0044] The present inventions are based on the hypothesis that spatial and temporal variations in T wave morphol- 
ogies and the relation between the depolarisation and repolarisation patterns will offer new measures of repolarisation 
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abnormalities. It was aimed to define a set of parameters that would (a) quantify such abnormalities, (b) be highly 
reproducible, (c) have sensitivity and specificity greater than the conventional measurements such as QTd in separating 
normal and clinically relevant abnormal electrocardiograms, and (d) be independent of problematic time-domain meas- 
urements such as the detection of the T wave offset 

5 [0045] ECG decomposition by several different methods is known. One example is taught by Acar B. and Koymen 
H in an article entitled °S VD-based on-line exercise ECG sig/ial orthogdhalizatibn",. (1999) IEEE Trans; Bibmed. Eng., 
46, pp. 31 1 1 -321 . The common approach is to transform the multiple lead ECG into another domain or to identify the 
dominant components of the recording. The analysis is subsequently carried out in the transformed domain -which may 
. have various advantages, e.g.. a high signal-to-noise ratio. ., 

io [0046] In the following description, the uppercase. bold letters are used for.matrices whereas the lowercase bold 
letters for vectors 

[0047] The analysis is based on the Singular Value Decomposition (SVD) of the standard 12-lead ECG It defines a 
minimum dimensional subspace that captures this ECG energy SVD was defined by Gdlub G H, and Van Loan C F 
(1996) in 'Matrix Computations', 3 rd edition, (The John Hopkins-University Press, Baltimore and London), pp. 70-71, as. 
75 [0048] If M is an 8-by-n matrix, with each row corresponding to a standard ECG channel (I, II, V1 , V2, V3, V4, V5, 
V6).and /? being the number of samples, then there exist orthogonal matrices 

U = [u„ ear- and V = [v, , .,v „] e 

20-}. j - . ' 

such that 

25 I^U^MV = &¥(<T J <T^)€«"* - f=triin{w,«) 

where a n £ o 2 £ > c p £ 0 

[0049] The columns of U are referred to as the left singular vectors, whereas the columns of V.are referred to as the 
30 right singular vectors, a, are the singular values. F.u rthermore, if a, > £ <s r > a r+1 =• ••= o p , = 0. then 

rank(M) = r, 

35 , s 

miiI(M) = span{v r4l , ^v,.} 
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where ranged) is the minimumdimehsionai space which captures the 1 2-iead ECG energy As shown in. the following 
section, no significant errors are. introduced by restricting the dimension, of this space to 3. The singular values are 
measures of how much ECG energy exists along the corresponding vector u 

[0050] The. signal representations in the minimum (rndimensional) decomposition sub-space and the corresponding 
left singular vectors arc used to derive the new parameters (Figure 1 ) In the rest of this text, M designates the ECG 
data matrix. Each column of M corresponds to a sampling instant and. each row corresponds to a different ECG channel 
Because of the algebraic interdependence eight of the standard 12 ECG channels are used, namely I, II, V1, V2, V3, 
V4, V5 and V6 The S VD .of M is performed as. 
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MV=U r MV 



I em M and is diagonal , U ett 



8xJ 
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wherer^" and is diagonal, and fen*" . It has been shown previously- that: 99% of the. ECG 'energy can be represented 
in a 3<fimensionalminimum subspace (see the paper mentioned above by Acarand Koymen, 1999)/Hence the* ef- 
fective rank of M is 3. This minimum subspace is spanned by the columns of ik*" Let Sbethe projection of Montou., 
s =f) r M Figure 2 shows an example of the input and decomposed (projected) signals. The transformation of the 3 
.dominant-decomposed signals back .into, the original ECG domain is equivalent to morphological.filtering of the ECG 
in its original domain 

. [0051] An approximate QRS and T wave detection is performed on the.3 , S2> s 3 ) decomposed signals that contain 
mostof the. energy; Let ; .'• r 



10 . 



= [s 2 *(',) *M] 

[0052] Figure 3 shows $ 3D and E 3D for a single beat The R-wave end point is assumed to be the first point after the 
maximum of E 3Z> where E 3b falls below 70% (arbitrary threshold) of its maximum value That instant is disclosed as 

20 t HE Although t RE \s nottheactual R-waveend point, it serves the purpose Similarly, the 70% point before themaximum 
is marked as t' RB The QRS complex is assumed to start 48 msec before T RE arid end 48 msec after ^(arbitrary 
limits) These two points are marked as f'^and t' RB respectively The T wave peak t TP is assumed to be the maximum 
point of E30 after i' RB The approximate T wave starting point t T$ is taken, to be 173 of (tip ifn£ after t'^. Figure 3 
shows the positions of t' RB V RB V RB t' RB t TS and t m Since this study is only concerned with single beats, there is ho 

25 need to choose an end point for the search of the T wave peak. However, such a point can easily be selected based; 
on the instantaneous heart rate, when necessary. The. detection of the T wave end point, t TB is based on the path of 
the tip. of Figure 4 shows the path that the tip of s^ott) follows in the plane spanned by u-, and i^, for t T& The 
rectangular area is divided into 100 equal rectangular cells (arbitrary number) Each cell is assigned a weight equal to 
the number of inner data points (time instants.) This is a measure of the time spent by the tip of in that cell Let D r 

30 be the weight of the 1 th cell The cells with zero weight are discarded and the other cells are ordered in respect of P r 
Figure 5 shows D, values for a single beat Assuming that there are * cells with nonzero weights, D, £ Cfe < • < D K 
[0053] A threshold th = mean(D t ) +\iX standard^deviationlD,) is used, where \i =3 (arbitrary constant). If D f Z th 
forL< i< Kior Y< i< O, the earliest time instant at which the tip of s 2D , enters one of the cells Y\o Q(L< i< K). is 
set to be the' approximate T wave end point, t TB 

35 [0054] Since the aim of the algorithm is to quantify the T wave shape between f TO and t m rather than to measure 
the V RE interval, the approximate and arbitrary nature of the f re instant definitions is fully acceptable (as shown further) 
[0055] If, using the algorithm described above, t m which should not occur, u. is increased in steps of 0.2 until t TB > 
tjp. Similarly, if D K < th, u. is decreased in steps of 0 2 until D K Z th. Such cases are rare 

[0056] This T wave end point detection scheme is based. on the concept that the interlead relations do not change 
40 in the absence of the ECG signal Each, point on the 14 u^-plane corresponds to a specific interlead relation defined by 
the vectors u 1 arid u 2 Hence, each cell in the u^-plane represents a group of similar interlead relations. When the. 
repolarisation pattern ends, the ECG signal remains confined to a small set of such relations 
[0057] The decomposed signal is subsequently normalised with the maximum energy set to t. 
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[00.58] QRS complex and the T wave extracted by this algorithm result in decomposed data matrices S r and S QHS> 
ss respectively From both signals, a DC vector is subtracted as follows. 
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s£=0.25x{s, o (^ 

5 [0059] From now on, 

s r €^ ? ^ W'=:in -t n ) and g;^ 3 ^ -i^}: . ; 

to • •■ . ... 

• will denote the decomposed, energy normalised and Decompensated T wave and QRS complex 
[0060] The T wave is reconstructed from $ v which is equivalent to a morphological filtering: 

15 M r = 0S T = 0tJ r M r: 

sThe reconstructed T wave, M r , is. once, again decomposed by SYD: 

20:, : 

25 

is diagonal. 5^3^ £ r ^9t'- > and is diagonal^ T ** M 

. [0061]. The subscript T indicates; that we are dealing with the T wave only, and the superscript w denotes a matrix: 
which is reconstructed (morphologically filtered) Note that*, has two columns whereas p has three columns This is 
because the 3 rd decomposed signal in T wave decomposition has been excluded (Figure 6) 
so [006.2] The spatial variation descriptors are determined as follows, ij; is an 8-by-2 matrix Its columns, are the two 
most significant left singular vectors o1m„ 



0 Z r 
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Each of its rows is the reconstruction vector of the corresponding standard ECG lead (note that the reconstruction from 
the 2D, most significant, subspace of the decomposition Space means multiplying the decomposed data matnx by u r ) 
Let z denote the reconstruction, vectors, 

40 

45 (note that the ECG energy along the two orthogonal dimensions of the decomposed space n TJt and a„ are proportional 
to the corresponding singular values or and g t £ To guarantee that we deal with the. morphologies rather than 
the energy differences, the decomposition space is rescaled to equalise the energies in both directions. 



so _ ~ ~ r ir ~ 



= .[W; W. y/ W n Wy 2 W K J W^ 4 W K5 W K J 
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Each w y represents the reconstruction coefficients of the T wave of the j* channel of the standard ECG 
[0063] The angle between different w ; vectors is calculated: 
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6. =Z( W| , W/ ) V/,y€{7,//^l,K2,K3,/4,K5^6} ( /^7 
5 ' e[0\i80 a ] 

. [0064] The smaller^ the closer the reconstruction vectors of the I th and the ^ .ECG channels. It was observed in ■ 
the data of this study that the T-wave morphology in. VI is generally different than that of other channels, irrespective 
to of any clinical background, mainly due to the position of the V1 efectrode Figures 7a and 7b demonstrate the difference 
between a normal ECG and a HCM case. The projection of the T wave loop onto reconstruction vectors gives the 
particular T wave as observed in the corresponding ECG lead. It is seen that reconstruction vectors for normal ECG 
are closely grouped (meaning similar morphology), whereas. they are dispersed (meaning different morphologies) for 
HCM patients Note that the reconstruction vector of V1 in normal ECG is far from the others This is the reason why 
'is it is preferred to exclude the V1 reconstruction vector from the calculations of the T wave Morphology Dispersion 
descriptor (i e , the average of angles between all pairs of reconstruction vectors), to be explained in greater detail 
below. 

[0065] The descriptor, T wave Morphology. Dispersion.(TMD), is defined.as. the rriean of all 0 7 y excluding V1 ; 

20 

mmv = ^£3, *U € \umv3,V4 9 vs % ve} J* j 

25 TMD ts a measure of the spatial T wave morphology variation 

[0066] Since the ascending and the descending parts of the T wave are known to correspond to different facets of 
the repolarization process, descriptors TMD pre and TMD posP which, are defined in the same way as TMD with the 
'ascending' part of the T wave (t rs <,t<t TP ) used for TtiD pre and the 'descending? part (fa £ t< t TE ) used for TMD post 
were also computed 

30 [0067] Both of the .QRS and the T wave represented by S QRS , and S r follow an. approximate loop in the column 
space of fj t The orientation of the T wave loop is determined by selecting the unit vector e T1 , with the maximum T 
wave energy, subsequently the unit vector e r?2 perpendicular to with the maximum energy, and finally the unit 
vector 0 r3 perpendicular to both e T1 andej^ 

[0068] The descriptor Total Cosine RJo_ T (TCRT) is defined as the average of the cosines of the angles between 
35 efA ^nd s Q ^f) (columns of S QR £j for all i w'rthin [t RS! t R £ This is a measure of the vector deviation between the 
depolarisation and the repolarisation waves 



TCRT = T -^ r - ^cosfzte,, ,s Qlts (i))). 
Ire { rs t-t*s 

[0069] Figures 8a and 8b show examples of the QRS and the T wave loop onentations in a normal ECG and a HCM 
45 case TCRT measures the deviation between these two loops It, in effect, measures the. difference between the prop- 
agation directions of depolarisation and repolarisation waves Negative TCRT values correspond to large differences 
in the loop orientation In Figures 8a and 8b, it can be seen that the orientation of QRS loop and T wave loop are close 
to each other in the case ol a normal patient, whereas they are far from each other in the case of an HCM (hypertrophic 
cardiomyopathy) patient. 

so [0070] The T wave representation S r , is de^normalised, removing.the effect of energy normalisation and its projection 
onto span(e TtV e T2 ) is calculated orienting the T loop along e ri The rectangular area encompassing the T loop is 
divided into n>>1 (in the present implementation n=4900, arbitrary constant) cells of equal size. The loop is closed 
with a straight line connecting the end points and re-sampied with equal sampling steps of the 2D space This re- 
sampling assures that there is at least one sample in every cell that the loop passes through The numbers of cells in 

55 the loop area and in the outer area are counted and the descriptors, Percentage of the Loop Area (PL) and Percentage 
of the Outer Area (PO) are calculated PL is the proportion between the number of cells inside the loop and the number 
of all cells. Similarly, PO is defined for the cells outside the. loop 
Note that PL + PO < 1 because there are cells occupied by the. loop itself 
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[0071] The descriptor termed Lead Dispersion (LD,), is calculated using the decomposed, energy normalised and 
DC-compensated T wave S T , as it was after the DC-compensation step. * T£D (/) £ span(u„ u 2 ) V /<; t T$ follows 
a path in the u, tvplane, the T loop The rectangular area containing that path is divided into 100 equal cells (arbitrary 
constant). LD, is defined as the number of different cells that the path involves It serves as a measure of the temporal 
variation of the interlead relations, during T wave. ■ . 

[0072] The descriptor, Lead Dispersion^ (LD 2 ), is defined similarly using [the denormalised T wave. It is the number 
of different cells that the T looppasses through, excluding the straight line that was added to close the loop. The basic 
difference between LD, and LD 2 is that Lb, is calculated using the energy.:normalised decomposed signals (maximum 
ECG energy = 1 ) while LD 2 is calculated using the original decomposed signals, 

[0073] : These descriptors provide further methods of quantifying abnormalities of ventricular depolansation and re- 
polarisation, and thus relate to further inventions disclosed herein 

[0074] The analytical system was implemented on a standard personal computer with Pentium 1 33MHz CPU and 
BOMB RAM using Matlab Version .5.2.0 (The MathWorks Inc., 1998). The system was tested with standard I2^lead 
ECGs recorded by the MAC VU Electrocardiograph (Marquette Medical Systems, Milwaukee, Wisconsin, USA) 10 
second recordings with 500Hz recording rate were acquired and the so-called median beat was obtained for each 
channel of the recording These median beats, sampled at 250MZ, were used in the analysis 
[0075] The inputs to system are 8, time-aligned, median beats, each one being a representative of the ECG mor- 
phology in the corresponding standard channel . The. main time consuming computations are the area calculations 
which involve a recursive algorithm, also implemented in Matlab, On the average, the computation of all parameters 
for a single recording takes 177 seconds If excluding the area related parameters (PL and PO), the analysis takes 30 
seconds per recording Matlab commercially available library without any modification is presently used A purpose 
built library would increase the performance of the system cpnsiderably. 
[0076] Three sets of. ECG recordings . were.iused in the study: 

(a) Standard resting 12-lead ECGs recorded in 1100 normal healthy subjects, 91 .3 male, aged 33±1 2 years, range 
10-81 years 

(b) 10 supine resting ECGs were recorded in each of 76 normal healthy subjects, 37 male, aged 38+10 years, 
range 13-59 years. In each individual, the serial ECGs were recorded one immediately after another using the 
same electrode attachments and without the subject moving during the whole recording session Data acqu.srt.on 
of each recording lasted 10 seconds and, including the electrocardiograph handling, each series of the 10. ECGs 

was obtained within 3 minutes . . 

(c) Using the same recording strategy, 10 supine resting ECGs were recorded in each of 63 patients with hyper- 
thropic cardiomyopathy (HCM), 44 male, aged 39± 14 years, range 12-71 years 

[0077] Using a research version of the commercial QT Guard software package. (Marquette Medical Systems, Mil- 
waukee Wisconsin, USA), several conventional descriptors of repolarisation. patterns were calculated for each ECG 
for comparison. This software aligns all beats with respect to the.Q wave onset and was programmed to use the 
downstope inflex tangent method to detect Ihe T wave offset The following conventional QT mterval and T wave 
parameters were considered: 

i Global QT Dispersion (G-QTd) = Max (QT interval in 12 leads) - Min (QT interval in 12 leads) 

ii Precordial QT Dispersion (P-QTd) = Max (QT interval in 6 precordial leads) - Min (QT inteivaL in 6 precordal 

Hi Area QT Dispersion (A-QTd) = A11 of the 12 leads are assumed to have the same T wave onsel and offset 
points- The areas under the T waves are calculated and the points atwhich they reach 90% of the corresponding 
total area are marked for each lead The dispersion, (maximum minimum) of these markers over 1.2 leads is cal- 
culsted 

iv Global j to Tpeak Dispersion (G-JTpd). = Max [J to Tpeak interval in i 2 leads) - Min (J to Tpeak interval in 1 2 
leads) 

v Corrected QT Interval (QTc interval). Bazett formula corrected maximum QT interval in all 12 leads 

Principal Component Analysis of 12 lead T waves is also incorporated in the QT Guard packages The B components 
with associated eigenvalues are obtained. Each eigenvalue is a measure, of the significance of the corresponding 
component If s,- denotes the. eigenvalue associated with the i* principal component, the following descriptors are 
calculated. 
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PCA ratio 1 (PGA,) - - 1 Jj=x 100 . 




*0 vii. 



VIII 



PGA ratio 2 (P0A g ) = j X 1 Q0 



S 3 

PCA ratio 3 (PCA 3 ) = — x 100 



Note that these parameters. can also be calculated by.SVD, assigning s ; = 

[0078] To find out whether the new methods assess. ECG qualities additional to.the conventional parameters, the 
correspondence between the new and the conventional parameters was investigated. AH new and conventional pa- 
25 rameters of 1100 ECGs acquired from normal subjects were used for this investigation and to calculate the Pearson 
Product-Moment correlation coefficient (Statistica Package, Release 5 1) between the new and the conventional pa- 
rameters and the ages of the subjects. 

[0079] The reproducibility of all the parameters was assessed based on the variation of the measurements between 
serial ECG recordings from the same individual of the populations of 76 normal and 63: HCM subjects The ratio of the 
30 individual range to the total range -was calculated for each patient and each parameter More precisely, for a fictional 
parameter X, this ratio fl*^ for normal subject j is ecjual to 



max(x*)-min(*;) 



40 where N=76 for ou r data set of 76 normal healthy subjects and X% is the value of descriptor Xin m lh ECG of the subject n 
[0080] The values fl£ were obtained in a similar way and for each descriptor. X } the means and standard deviations 
of values flj^ f and ftj^j were calculated and used tb compare the reproducibility of aill the descriptors 
[0081]. Weafso investigated the univariate and multivariate distinction between normals and HCM, using all of the 
new and conventional descriptors The parameters were compared on. the basis of the significance in discriminating 

45 these two groups and in terms of specificity and ! sensitivity The populations of 76 normal aind 63 HCM subjects were 
used and for each descriptor the mean values of 10 recordings were considered 

[0082] Individual parameters were firstly used in a univariate analysis The normal and HCM groups were compared 
using non^parametric ManrVWhitney test implemented, using an. in-house, written software according. to the original 
description (see Maim H B. and Whitney D R. (.1 947). 'On a test of whether one or two random variables is stochastically 
so larger than the other", Ann Math. Statistics, 18, pp 50-60 for more information in this regard) .P-value < 0 05 was 
considered as statistically significant 

[0083] The Receiver Operator Characteristic (ROC) curves which show the dependency of specificity on sensitivity were 
c;aiculated1or each individual parameter usingan iri-house software package (see Hnatkova K., Polomeeki-J.D: -Gamm A 
J , Malik M. (1.994). "Computation of multifactorial receiver operator and predictive accuracy characteristics", Cornp, Meth 
55 Prog Biomed., 42, pp. 147-156, for more information in this regard). The area under the ROC curve (reported as a per- 
centage) was used to characterise the predictive power of eaich parameter Independent of fixed sensitivity Jevels. 
[0084] Multiple regression analysis was used to assess the relative performance of individual descriptors in discrim- 
inating the HCM from normal subjects. The dichotomy limit of each parameter was set to the mean of the average 
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values of the norma! and the HCM groups. Multiple regression models of different orders were calculated (Statistics 
package, Release 5. 1) in a backward stepwise manner, by excluding the least significant variable at each step until 
the p-values of all surviving parameters were below 0 05 

[0085] The most significant parameters, which were identified by the multiple regression analysis, were used for the 
.calculation of. multivariate ROC. curves, The combinations of two and four parameters were used with the decision 
rulesof atle^st 1 of 2 positive" a ; hd at least 2 of 4 positive; respective^." ": 
[0086] ' The. resultsobtained using the new and conventional descriptors were as follows: 

[0087] Table 1 gives the Pearson Product-Moment correlation coefficients between the- conventional and the new. 
descriptors and the age of the patients!. None of the new or;conv^ 

the age ( I r I <<© 1 6 for all parameters): the absolute valuaof all thecorre lation coefficients between new descriptors 
were < O.S except for: TMD/TMD^: 0.9.1, TMD/TMD pre : 0.93, TMD pos /TMD pre : 0 79,- PUPO. -0 94, PL/LD 2 : -0 54, 
PO/LD 2 . 0 .50.. . 

[0088] The absolute values of the correlation coefficients between the conventional and the new parameters, were 
all < 0 5, except: TMD/PCAg: 0 552, LD2^CA 2 . -0.562 

[0089] The reproducibility of the conventional descriptors was generally poorer than that of the new ones (Figure 9) 
with the exception of PC^ and PCA 2 , that had reproducibilities .similar to those of the new descriptors 
[0090] Table 2 shows the mean values of .e # observed in 1100 normal subjects and confirms the reasons for excluding 
V1 lead from spatial variation descriptors An increased spatial variation of T wave morphology, in HCM subjects was 
observed in all of the three spatial variation descriptors, i.e TMD; TMD pre andTMD post . The meahTCRT was negative 
for HCM subjects arid positive for normals. The PL was larger (or normal subjects whereas the PO was smaller The 
mean value of LD, (as well as of LD 2 ) in norinaj and HCM. groups were close to each other 

[0091 ] Table 3. shows the comparisons between the descriptors in normal and. HCM subjects While all of. the de- 
scriptors strongly differentiate t?etween both groups, sorne descriptors have substantially lower p-values than others 
The QTc interval. and P-QTd offer the. most sign if icant. univariate differentiation among the conventional descriptors 
. However^TCRT, TMD and TMb^outperlormed all pf the conventional descriptors, while TMD pre had a p-yalue close 
to that of QTc interval, wtiich is the best ainong the conventional descriptors. 

[P09.2] Table 3 also shows the area under the univariate ROC curve for each descriptor. The results confirm the 
statistical comparisons: TCRT TMD and TMD post have areas above 90%, QTc interval has the largest area (85.6%) 
among the conventional descriptors 

[0093] Table 4 shows the p-values of each parameter in a succession of multiple regression models pf different 
orders The descriptors TCRT, TMD pre . P-QTd and QTc interval survived throughout the successive multipleregression 
models performed in a backward f stepwise fashion When excluding TMD pre and-TMD post , the final significant param- 
eters surviving the backward stepwise multiple regression analysis were TCRT, TMD.. P-QTd and QTc interval, with, p- 
yalues of 5.75X10- 8 , 0 011, 5 93x10^, $.61*10* respectively In both cases, TCRT outperformed all of the other 
parameters in all orders of multiple regression analysis 

Table 1 



Pearson Product-Moment correlation coefficients between the new descriptors, the conventional descriptors and 
the age of the subjects. • , . 



CORRELATION COEFFICIENTS 





TMD 


TMDp^ 


TMD pre 


TCRT 


PL 


PO 


LD, 


LD 2 


TMD 


1 00 
















TMD po s t 


0.91 


1.00 . . 
















0 93 


0 79 


1 .00 












TCRT 


-0 01 


I 0 04 


0.05 


1 00 










PL 


-o iO ' 


16 


: -0 18" 


^008 


1 00' 








PO 


0.O9 


.0 14 


017 


0 08 


-0 94 


1 00 






LD, 


0 08 


0 05 


0.07 


.0 00 


0 14 


-017 


1 00 




LD* 


:0 .30. 


. -o.i:a 


:0 17 


0 20 


-0 54 


050 


-012 


1.00 


G-QTd 


0 08 


0 10 


0 01 


O 03 


O.04 


-002 


-0 01 


-0 08 


P-QTd 


016 


0 16 


0 09 


0 00 


0 03 


-0 02 


.0 0.0 


-0 15 


A-QTd 


.0 10 


011 


0.00 


-0 08 


004 


-0 03 


0 01 


-0 13 
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"*'-■' Table T (continued) 



CORRELATION COEFFICIENTS 




TMU 


T MDp OSt 






PL 


PO 




LDo 

2 


G-JTpd 




o. iy 


n no 

• u.uy. 


«n 1 9 




-0^22 


.0.02 


. -0.48 


PCA! 


0.25 . 


0.1.7 . 


: :.0:22 


0.09 ,' 


0.06 


-0,05 


0.15 


-0.24 


PCA 2 


. 0.55 . 


: 0.46 


. .0.46 • 


.. -0.14 . . 


: . 0.04 . 


-6.06. . 


0.05;. 


; -0.56 


pca s : 


.0.16 ' 


.6.15. 


0.14 


• o oo 

0.03 


-U.ll 




n 11 


-0.12 


QTc 
interval 


o;o7 


0 07 


0 08 


-0,05 


-0 17 


016 


0 04 


0 04 


AGE 


-0 06 


-0.11 . 


:• 0-00 


.0:09 


. -0.06 


0 01 


0:05 


0.07 



AVERAGE d ti VALUES 



Channel 


I 


u 


VI • 


V2 


V3 


V4 


V5 


V6 


I 


0000° 


6799°+ 
7.997° 




22 47°+ 
24!82° 


10 65°+ 
10.79° 


4,648°± 
6.553° 


4 451°± 
4 603° 


5.966°± 
6:218 


II 




0(KK)° 




27 44°± 
. 25.79° 


15.. 43°+ 

12:79° 


7.575°+ • 
9.010° 


4 70l°± 
6 317° . 


4.180°+ 
6:015° 


VI 






o:o6o° 


IPIil 




»^ 

i mm 






V2 








0600° 


1.2 94°± 
21.72" 


.20 73°± 
24.39" 


24 86°+ 
24.97" 


27.65°± 
25.42° 


V3. 










0 000° 


8 6K5"± 
.8 931° . 


12 84°± 
10.70° 


15 64°± 
12.31° 


V4 












0000° 


.4 406°+ 
5.460° 


7 209°± 
8.129° 


V5 














0000° 


2.939°± 
3.893° 


V6 
















0.000° 



Tabic 2 Average © aj values in the set of 1 100 normal subjects The shaded boxes contain the 8i.vi and 
Bvi., values that are greaiter than the others 



Table 3 



The mean values and the standard deviations of all parameters, the Mann-Whitney Test results and the area under 
univariate ROC curves of the separation between normal and. HCM subjects. 


SEPARATION BETWEEN NORMAL AND ABNORMAL ECGs 


Parameter 


Normal n=76 


HCM h=63 


Mann- Whitney Test 


Area under ROC 
Curves 


Mean.± SD 


Mean ± SD 


P-Value 


TMD 


10 72 ±4.784 


41 10 ±26.85 


! 2 818X10* 18 


90 1% 


TMD post 


6 141 ±4.462 


36. 68 ±27 49 


2 289X1 0" 19 


91 1% 


TMD pre 


8. 682 ±4: 585 


4214±32 62 


1 605X1 0-1 3 


85.1% 


TCRT 


0 522 ±0 274 


-0 351 ± 0 522 


3.548X10" 19 


909% 
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Table 3 (continued) 



SEPARATION BETWEEN NORMAL AND ABNORMAL ECGs 




Parameter 


Normal n=7S 


HCM n=63 


Mann-Whitney Test 


Afpo under ROC 

. , Curves 




Mean ± SD * 


i Meani.ou 


i value 




PL ■ 


0.671 ±0.085 • 


c\ cno a.- r\ i AO ' ' 
U.oUo±U. ; 14Z 




64.3% , 


. PO 


. .0.273 ± 0.072 - 


0.32o±y.ll5 


•' fi^l V 1f)3. 

O.UC> 1 A IU • . 


• 65>/o 


LD n 


36.40 ±1.163 


34.81 ± 3. 15/ .. . 




71.8% . 


LD 2 


724:5 ± 346.1 


604.9±4ob.i 


D. / O / , ^ I U 


67 4% 


G-QTd 


19 97 ±11.62 


36 55 ± 1 o Bo 




77.5% 


P-QTd 


10 79 ± 8 776 


0"7 fl7 _1_ -i O CO 

27 87 ± lo by 


D D II A IU 


80.6% 


A-QTd 


13 70 ±8.564 


24. 38 ± 12 23 




76 8% 


G-JTpd 


32 53 ± 12,18 


45 9b x ZV F o\ 


aZ. 400A 1 W 


70.8% . 




680 0 ± 226 3 


481 4±245.8 


6 698X10" 8 


76 7% 


PCA 2 


15. 56 ±5 162 


23 56±1Q85 


. 9.886xl0f ' 


. 74 4% 


: PCA 3 


4.826 ±2,373. * 


7 765 ± 4. 235 
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Table 4 Significance levels of the parameters in different orders of multiple regression models, 
calculated in a backward stepwise fashion by excluding the least significant parameter at each 
step 
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[0094] The multivariate ROC curve involving the descriptors TCRT. TMD pre , P-QTd and QTc interval, had an area 
j- of 98 4%. Using the conventional and new descriptors separately in bi-varjate' ROC curves, we obtained areas of 95.5% 

\ for TCRT and TMD pre , and of 91 6% for P-QTd and QTc interval 

I [0095] Hence, TCRT, TMD pre » P-QTd and QTc interval are mutually independent separators of normal and HCM 

5 EGGs of which TCRT- is by far the strongest. 

[0d96] The neW parameters proposed in this specification are defined using the decomposition space and aimed at 
. the description of the temporal and spatial variations of ventricular repolarisation. The descriptors TMD, TMD pfe and 
TMDpos, reflect the interlead : morphblogical variationsof the T wave patterns* that is the spatial variations. The area 
: . related descriptors, that is PL, PO, LD, and LD 2 , characterise the temporal variations, TCRT introduces the concept 

io of; comparing the global, wavefront directions of the depolarisation and repolarisation processes. 

[0097] The original hypothesis was' that, compared to norma! ECGs » the s P atia! ^ temporal variation of T wave 
morphology are increased and the depolarisation. and the repolarisation vectors are more different in pathological 
recordings, such as in HCM patients. The statistical comparisons of this study verify this hypothesis. 
[0098] The mean of PL is higher and PO is lower in normal than in HCM subjects. This suggests that the T loop is 

75 relatively smooth and connected (not crossing Itself) in normal ECGs than in HCM EGGs On the other hand, despite 
the significant difference, the lead dispersion parameters (LD, and LD 2 ) have similar mean : values for both the normal 
and HCM subjects This suggests that the loop itself is not discriminative The loop lengths were similar in both groups 
The discrimination by PLand.PO seems to be due. to the disconnected and narrow loop (inner area is similar to a strip), 
rather than due to an increased irregularity of its shape 

20 [0099] The cfiange of sign of TCRT between the normal and the HCM subjects provides a clear distinction between 
the two groups (the negativity of TCRT shows an increased deviation) This is in agreement with the original hypothesis.. 
The repolarisation and depolarisation waves do differ in terms of their principal direction in a 3-dimensional time-or- 
thogohal space. Since the mean difference between normal and abnormal ECGs is52 0 , the descriptor does not merely 
reflect T wave inversion that would result in the difference near to 1 80° 

25 [0100] TCRT, TMDp re , P-QTd* and 'QTc interval were the only parameters that survived throughout the backward 
stepwise multiple regression analysis comparing normals and HCM subjects and TCRT was the best throughout the 
test This verifies that the new spatial variation parameters and TCRT are. very potent descriptors of repolarisation 
abnormalities 

[0101] In this preferred embodiment, aJI of the new descriptors are defined using the decomposition space This 
30 provides an inherent immunity to noise and avoids the inaccuracies associated with time domain measurements, that 
are common in QT interval related descriptors that depend on T wave offset determination The independence of time 
domain measurements makes the new descriptors highly reproducible, which is very important for their potential clinical 
applicability Among the. conventional parameters, only PG ^ and PCyV have a reproducibility in the same order which 
is again due to avoiding the. time domain measurements. 
35 [01 02] The weak correlation between the new and the conventional parameters shows that the new concepts quantify 
different properties of the ventricular repolarisation Furthermore, the two new concepts, that. is the spatial and temporal 
variation can be identified by strong correlations within each group, and by weak correlations across the groups The 
concept of TCRT is different from both spatial and temporal variations and the descriptor does not correlate strongly 
with any other descriptor. 

40 [01.03] The relatively poorer reproducibility of the T loop area related parameters is due to the algorithmic problems 
An open loop may result from baseline wander, as well as ST-segment elevation/depression A straight line was used 
to connect the ends of the loop, which is not necessanly the best approach An alternative may be to connect the ends 
of the loop and its centre of gravity or to transform the a, c^-plane creating a closed loop. It is alsd possible that the 
loop crosses itself, resulting in more than one surrounded area. In the study, the inner area was defined as the closed 

45 area neighbouring the beginning 6f the loop, ignoring the "pockets' It is an open question whether the existence and/ 
or the area of these 'pockets' is of any significance. The poor performance, of the loop related descriptors in differen- 
tiating normal and abnormal ECGs may well be due to these problems 

[0104] The arbitrary choice Of constants, used in dividing the plane of ECG into equal size cells, have an infFuence 
on LD V LD 2 , PL and PO calculations as well as in approximate T wave offset detection They define the precision of 
so these descriptors. Increasing these. constants would increase the precision at the cost of increased computation time 
However, the precision is also restncted by the ECG sampling rate which determines the smallest distance between 
two consecutive ECG vectors. Unreasonably decreased cell size (increased constants) would also degrade the per- 
formance of T wgve offset 'detection 

[0105] On the other hand, the T wave onset/offset definitions may have an influence on the temporal variation de- 
55 scnptors but do not affect the others Setting the constant [iinT wave offset detection to 3 is an appropriate choice 
The algorithm readjusted that value in 91 of 11.00 normal ECGs The QRS onset/offset definitions , on the other hand, 
are robust and able to handle wide QRS complexes. However, the choice of 70% threshold in determining the region 
of QRS used in TCRT calculation is important A too low threshold may result in a too general estimation of the QRS 
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loop orientation, whereas a high threshold may/misinterpret the .orientation of the-depolarisation wavef ront vector. 
[01 06] Only the principal direction of the EGG vector during X wave, was used in.TCRT calculation, that is the 
other lower energy components, e T? and b tz were ignored The average ratio of the energy along the second com- 
ponenttothat of the first was 0. Ufornormals and 0.22i6r HCM patients. Thisshows that the T loop generally resembles 
a narrow ellipsoid and it is the direction of this loop that is of interest. Using e T2 and/or e Tj3 would not improve the 
concept of TORT rnaihiy due to a decreased iio'ise immunity. There is no ambiguity cri the e' Ti1 definition because the 
Decompensation ensures that.e T1 hastha correct sign. 

• [0107] In conclusion, therefore, -.it can be seen that the new descriptors of repolarisation patterns descnbed. in this 
specification have several important qualities. 

• All of the new descriptors can be assessed in a minimum dimensional space constructed by SVD of 1 2-lead ECG 
This provides a built-in immunity to noise. 

• None of the new descriptors require accurate time domain interval measurements This makes them, more repro- 
ducible than the conventional QT interval related descriptors 

is . The new descriptors assess different ECG qualities than the conventional parameters 

• The spatial variation and wavefront direction, descriptors can discriminate between normal, and abnormal EGGs 
... substantially better than the conventional descriptors The wavefront direction descriptor (TCRT) is by far the 

strongest of all considered in this study 

so ,[0108] In a second case study, the effect of changes in the autonomic tone of a patient was studied using known 
descriptors and some of the new descriptors described in the first case study. The findings were as follows: 
[0109] In the early thirties, Wilson et al (Wilson FN , Macleod A G , Barker P S,, Johnston F.D , The. determination 

• and the significance of the areas of the vontricular.deflections of the electrocardiogram A. Heart J 1934; 10:46-61, 
Wilson F.N , Macleod AG., BakerRS* The T Deflection of the Electrocardogram. TR A Am Physicians 1931;.46.29) 

26 mtroduced the concept of algebraic sum of the areas under the ventricular deflections of the electrocardiogram (the 
net QRST area) The resulting vector quantity, called 'the ventricular gradient (VG) was proposed as independent of 
the sequence of ventricular activation as long as the ventricular recovery properties remained constant Therefore it 
was believed that the VG could help distinguish T wave changes following changes in the activation pattern (isecondary' 
T wave changes) from those due to myocardial damage ('primary' T wave changes) 

30. [0110] The original idea did not evolve into a clinically useful tool both because of technical difficulties with the meas- 
urement, and because of data challenging the quantitative independence of the activation sequence.. As discussed in 
relation to the first case study, a new descriptor of the wavefront direction of repolarisation has been proposed It 
quantifies the difference in. the global direction of the wavef rbnts of. the . depolansation and of the repolarisation as an 
average of the cosines of angles between main depolarisatiqn and repolarisation vectors in a minimum dimensional 

35 subspace derived! rom theindependent leads of the electrocardiogram (total cosine R to T, TCRT) Although measured 
in the optimised 3D space that contains the maximum :energy of all. ECG leads, this new descriptor advances the 
classical concept of VG 

[0111] Ventricular gradient is a vector which gives the direction and magnitude of the electrical forces produced by 
a lack of uniformity in the .duration of the excited state, it points from the region in which the average length of systole 

40 is greatest toward the region in which it is least. Once considered a 'fundamental quantity in electrocardiograph/ the 
VG was gradually forgotten due to both uncertainties about the validity of its. concept as well as to technical difficulties 
with its manual calculation from QRS.T time integrals before the personal computer era Today the VG is hardly men- 
tioned in modem textbooks of electrocardiology with comments such as that 'the most exciting thing about the ven- 
tricular gradient is its name' The. TCRT has been demonstrated to be more reproducible and to separate normal from 

45 abnormal ECGs better than several repolarisation parameters including the dispersion of the QT-end and JT-peak 
intervals and the corrected QTc interval. TCRT has also been found to contain independent predictor value of mortality 
and arrhythmic complications after myocardial infarction. 

[0112] These observations suggest that the concept: of VG might, have been neglected prematurely. Havmg-this in 
mind, several studies were initiated researching the prqperties of VG and of its modem and more precise counterpart 

so m the present study, the effects of basic autonomic provocations on VG were investigated The study assessed the 
effect of postural changes and autonomic provocative manoeuvres on the direction and magnitude of the spatial VG 
and on TCRT in healthy subjects. We compared the effects of postural and autonomic provocations on VG descnptors 
With the effect on conventional repolarisation parameters, namely the QT interval duration and QT dispersion 
[0113] The study population consisted of 40 healthy subjects, 31 male, median age 33 years, mean age 33.1+7 3 

55 years range 18-56 years, with no history of heart disease and with normal resting 12-lead ECG None of the subjects 
was takingany autonomically active medication and before the test, the subjects were instructed to refrain from smoking 
and from caffeine intake 
[01 14] The procedure was as follows. 
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[0115] Following 10 minutes supine rest in a comforting temperatured and dimmed room with a. low level of back- 
ground noise, the subjects performed the following tests. 

[0116] Postural changes (32 subjects), resting supine position, followed by sitting, unsupported standing, supine, 
and standing position, 4 minutes in each position (total of 20 minutes), with abrupt transition between the separate 
positions.. . . . : ' ' - .... 

[0117] Valsalva manoeuvre (30 subjects): forced expiration into the mouthpiece of mercury manometer maintaining 
a constant pressure of 40 to 50 mm Kg for up to 1 minute. Each manoeuvre was performed a times in supine position 
with 4 min periods between the tests, and 3 times in unsupported standing position ateo with 4 minutes periods between 
. the manoeuvres. 

101 1 8] Sustained Handgrip (8 subjects): the maximum force of contraction was determined in each individual follow- 
ing which each subject maintained 30% of maximum forcefor 5 min. Each test was performed twice in supine position 
with 4 minute periods between the tests and twice in unsupported standing position again with 4 minutes of rest between 
the tests. 

[0119] The data were recorded using continuous 12-lead digital ECG (250 Hz sampling rate, 1 2 bit A/D conversion) 
for the complete procedural section for each subject without any Ibss of signal using a digital recorder (SEER MC, 
Marquette Medical Systems, Milwaukee, Wisconsin, USA) ECG data were storeci.in separate 10 sec portions The 
individual tests were organised in a synchrony with the recorder in order to identify each 10 sec ECG. sample within a 
specific phase of each test 

[0120] From each lead of each 10 sec ECG sample, the so-called 'median beat* was constructed representing the 
ideal GRST complex of the given ECG. Data analysis of ECG patterns was based on these median complexes 
[0121] The magnitude and angle of spatial VG were calculated in the following ways From the median ECG beat, 
the area of the QRS complex and of the T wave in each of the 12 leads were calculated using the ECG Research 
. Workstation Software Package Version 1.0, by.Marquette G.E. (Milwaukee, Wl, USA). The QRS and T wave areas in 
the orthogonal X, Y, 2 leads were derived from the areas in the 8 independent leads (I, II, V-,, V 2 > V 3 , V4, Vg and V 6 ) 
using the inverse Dower matrix. The magnitude of the spatial VG was calculated as 

;VG m ^ ((QRS X + TJ 2 + (QRS y + T y f + (QRS Z + iff 

where QRS,, and T^ are. the areas of the QRS complex and of the T wave in the- orthogonal lead W, respectively The 
angle VGa of the spatial VG was calculated as the spatial angle between vectors. originating in the centre of the co- 
ordinates with the final points of [QRS x ,QRS y , QRSJ and [T^Ty/TJ, respectively. 

[01 22] As described in detail in the description of the first case study, the total cosine R_ to^T (TCRT) was calculated 
The median beat was represented in a minimum dimensional subspace using singular value decomposition of the 
standard 12-lead ECG The TCRT was defined as the average of the cosines of. the angles between the main QRS 
and T vectors in the 3-dimensjonat reconstructed subspace. In effect, TCRT measures the difference between the 
propagation directions of the depolansationand repolarisatibn waves, with smaller (and negative) values corresponding 
to greater difference between the two wavefront directions. 

[01 23] The median beats of all ECGs were analysed automatically using the QT Guard software package (Marquette 
G.E ) A common onset of the Q-wave in all. leads was identified and the offset of the T wave in each lead was determined 
by the downslope inflex tangent method For the purpose of this study the maximum QT interval, the global QT dis- 
persion (QTd', maximum Qt interval in 12 leads - minimum QT interval in 12 leads) and the RR interval were taken 
from the results provided by the QT Guard package 

[01 24] The first 2 minutes of 1he supine rest recordings were ignored in order to achieve fully stabilised steady state 
The recordings obtained during the final Sminutes of 1he 10 minutes of supine rest were used toderive baseline values 
for each parameter and to investigate their mutual correlation as well as correlation with heart rate To investigate the 
correlation, averaged supihe resting values of individual subjects were considered. 

[01 25] The mean vaiues of each parameter for the separate, positions and autonomic manoeuvres were calculated 
and compared by paired t-test and one-way within subjects, (repeated measures) analysis of variance (ANOVA) with, 
post hoc comparisons using Scheffe test (Statistica, Version 4.00). All values are expressed as mean ± standard error 
of the mean (SEM) Statistical significance was defined as p<0 05 

[0126] The results of the study were as follows. The correlation coefficients between VGa, VG, TCRT and RR. interval 
during steady-state supih e "position are shown in Table 5. There was a significant, correlation, between the angle and 
the magnitude of the VG, and between, the angle of the VG and TCRT While both VGa and VGm were significantly 
correlated to the RR interval, there was no significant correlation between TCRT and RR (Table 5) 
[0127] Postural changes significantly decreased both VG and TCRT While VGm and TCRT Were significantly de- 
creased in sitting and further in standing position position, VGa was increased in sitting and was further increased in. 
standing position (Table 6). 
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[01 28] As expected, the RR interval was significantly shortened in sitting and further in standing position. The max- 
imum QT interval followed the changes of the RR interval (see Figures 10a and 10b), 

[01 29J VGa,VGm and TCRT were abruptly changed with transition from supine to sitting and from sitting to standing 
position. There was a general tendency for TCRT to be changed more abruptly when assuming each new position, 
compared to both VGa and VGm. 

[0130] VG m and TCRT were, significant ly decreased and VGa was significantly increased during the strain phase of 
Valsalva manoeuvre compared to preceding resting period both in supine and in standing position 
[0131] The RR interval were significantly decreased during the strain phase of Valsalva both supine and standing. 
QT max was sligthly but statistically significantlyshortehed during Valsalva in standing. (364±4 vs 371 ±5 ms, p=0.02) 
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position (Table 6). 

[01 32] QT dispersion was not changed significantly during Valsalva. 

[01 33] VGm and TCRT were slightly but significantly increased by handgrip in supine position compared to the pre 
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ceding resting period (57.8*4.1 vs 56 1±3.9 mVms, p<0.0001 for VGm and 0.63*0.06 vs 0.6,1±0.06. p=0.0007 for 
TCRT) Both descriptors were not changed significantly by handgrip in standing position. VGa was slightly but statis- 
tically significantly reduced by handgrip in supine (40 0±3.1 vs 40 8±3 1, p=0 01 ) and increased in standing position 
compared to preceding resting period (52.8±3 8 vs 51 .8+3.7. p=0.027) 
5 [01 34] Neither the RR interval, nor QTd were changed significantly by handgrip compared to preceding resting period 
in supine as Well as in standing position:' However,- QTitiax was slightly;. but .signifipahtly,shp.rtened by,handgrip.in 

standing position 370±5 vs 374±5, p=0 06). • • ■ ■/ ". - , 

[0135] The aim of this study was. to examine the effect of. established,autonomie tests.on the spatial ventricular 
• gradient and on a new descriptor of the wave direction of depolarisatbn and. repolarization; TCRT . 
io . [0136] The main finding is that both the magnitude and the angle of the spatial ventricular gradient, as well as the 
new wave direction descriptor, TCRT, react sensitively and rapidly to changes induced by postural and Valsalva ma- 
noeuvres As expected from their mathematical relatipn (i.e the algebraic sum of two vectors decreases as the angle 
between them increases) VGm followed, the opposite trend to VGa and was significantly reduced in sitting compared 
to supine and was further reduced in standing position 
is [01371 Even more impressive than the absolute magnitude was the speed of the change of VG and TCRT during 
postural changes and Valsalva manoeuvre Statistically significant changes in TCRT were detectable already in the 
first or second 10-sec median ECG beat after standing up from supine or sitting position, or lying down from standing 
- position (Figures 10a and 10b) . Considering the fact that during the first 10-sec recording in each new position, the 
• . actual postural change took place (taking approximately 3-4 seconds), !! is apparent that TCRT is an extremely rapidly 
so responding parameter TCRT appears to be able to respond more rapidly to autonomic modulations than both VGa 

rai38] 3m previous studies have described significant reduction of VG m in standing position.compared to supine, as 
Well as during the strain phase of Valsalva manoeuvre. The effect of all phasesof Valsalva manoeuvre on repolansation. 
descriptors; however, can be assessed precisely, only on beat-torbeat basis analysis. 

25 [0139] It is difficult to comment on the differencesin the effect of Valsalva manoeuvre andsustamed handgrip on the 
wavefront direction descriptors Althoughhandgrip in supine position statistically significantly increased VGm (56.1+.3. 9 
vs 57 8±4 1 mVms, p<0 0001) and TCRT (0.51±0.06 vs (1 .63+0.006, p=0,0007) and decreased VGa (40. 8±3.1 vs 
40 0+3,1 , p=0.01 ) the differences were much smaller than those induced by Valsalva manoeuvre.and postural changes 
The handqrip test however, is known to be of limited sensitivity and specificity 

so [0140] We found statistically significant correlation (r = -0 7.8, p < 0 0001) between VGa and TCRT during steady- 
state supine conditions The spatial VGa and TCRT appear to quantify the same physiological phenomenon, namely 
the difference in. the spatial direction of the. wavef ronts of the depolansation and the repolansation, in a three<iimen- 

[oT4n $P VGa was positively correlated with heart fate (r = -0 39, p=0.026 for the VGa/RR relation). Ontheother hand, 

35 although TCRT and the heart period clearly followed the same trehd during postural. changes and.Valsalva manoeuvre 
(see Figures 10 and 11), there was- no correlation between the two parameters during steady state supine positioner 
= 0 09 p = 0 61) This suggests that rather than TCRT being driven by the heart rate, both, parameters are under the 
control of a common factor operating during autonomic perturbat.ons but not during steady state resting conditions 
[0142] Previous studies have found significant correlation between VG and the heart rate On the basis.of this relation 

40 VG found limited use as a sensor for rate-adaptive pacing The gradient is increased by an increase in pacing rate 
and is decreased by exercise or other stress during fixed rate, ventricular pacing. Thus, at least in theory, the gradien 
cart be implemented in a closed loop rate-adaptive f eedback system, the gradient js decreased by physical or emotional 
stress, which leads to increase o! the pacing rate (at a controllable speed) and the latter increases the gradient to the 
resting level. The gradient was measured by the time integral (the area) underthe evoked R wave Important limitations 

« of the VG measured in this way prevented the widespread use of VG as a rate-responsive sensor: although it is a 
rapidly responding sensor, the maximum pacing rate is reached very early during exercise, so the proportionally of 
the rate-response is very moderate The ventricular gradient is alsodirectly rented to the ventricular mass and thick- 
ness therefore changes in the ventricular geometry in upright posture may lead to paradoxical changes in heart rate 
[0143] The data from this study suggest that the role of VG and especially of TORT as rate-responsive sensors in 

so implantable antiarrhythmic devices, and possibly also for automatic detection of potentially arrhythmogenic autonomic 
modulationsare potential applicationsof these descriptors. To achieve this. in a commercial product, it may be necessary 
to make certain modifications such as simplyfying the measurement algorithm of TCRT and to calculate, it from intrac- 
ardiac leads instead of standard surface leads . 

[0144] In some recent studies it was found that the beat-to-beat variability of the VG was significantly increased by 
ss myocardiaJ ischemia (see, for example Horinaka S, Yamamcto H, Tabuchi T, Takada M, Akabane T, Onoda M Yag, 
Ventricular gradient variability. New ECG method for detection of ischemic heart disease J Electrocardiol 1995; 2B 
177-183 and Horinaka S, Yamamoto H Enhancement of ventricular gradient variability during acute myocardial 
ischemia Int J Cardiol 1998, 65.-173-180) In another recent study it was found that an abnormal T wave axis in either 
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the frontal or horizontal plane was a. strong and independent predictor, of fatal and nonfatal cardiac -events in the 
general population older than 55 years: TCRT is able, to observe-such an abnormal T-wave axis. 
[0145] From the resufts of this study, it was found that the maximum QT interval followed the changes of heart rate 
during postural changes and was significantly decreased in sitting compared to. supine position and was further de- 
creased in standing position. Although statistically significant, the changes of several milliseconds of QT max in standing 
positionby Valsalva and haridgrip.are hardly of any clinical significance. . 

[0146] Although QT dispersion was significantly decreased in sitting and further in standing position compared to 
supine; the values are largely overlapping. It has been shown in the past that QT dispersion was significantly increased 
in standing compared to supine position in patients with syndrome X, while in one study a significant effect of posture 
on QT dispersion in healthy subjects was hot found The low reproducibility of both automatic and manual measurement 
of the QT dispersion is well documented. 

[0147] The role of the autonomic nervous system in ventricular arrhythmogenesis has recently been heavily inves- 
tigated (sec, for example, La Rovere.MT, Bigger Jr JT : Marcus Fl> Mortara A, Schwartz PJ Baroreflex sensitivity and 
heart rate variability in prediction of total cardiac mortality after myocardial inlarction Lancet 1 988, 351.478-484, Hohn- 
loser SH, Kiingenheben T, Loo A, Hablawetz £, Just H, Schwartz PJ . Reflex versus tonic vagal activity as a prognostic 
parameter in patients with sustained ventricular tachycardia or ventricular fibrillation Circulation 1994, 89.1068:1073, 
and Schmidt G, Malik M, Barthel P, Schneider R, Ulm K, Rofnitzky L, Gamm AJ, Bigger Jr JT, Sehomig A. Hear-rate 
turbulence after ventricular premature beats as a predictor of mortality after acute myocardial infarction. Lancet 1999; 
353:i 390-1 396). The autonomic modulationsof the dynamicity of ventricular repolarisation have been shown to contain 
diagnostic and prognostic information (see; for. example, Maison-Blanche P.Coumel P. Changes in. Repolarization 
Dynamicity and the. Assessment of the Arrhythmic Risk PACE 1 997, 20[Pt ll],261 4-2624, Nollo G, Speramza G, Grasso 
R, Bonamini R, Man.giardi L, Antolmi R: Spontaneous beat-to-beat variability of. the ventricular repolarisation duration 
J Electrocardiol 1992; 25.9-1 7; Sandrone G, Torzillo D, .Fundar6 C, Porta A, Danna P, Polese A, Maliiani A, Lornbardi 
F Spectral Analysis of RR and R-t Variabilities in Patients With Coronary Artery Disease; A N E. T998; .3(3):237-243; 
Gang Yi, Xiab-Hua Guo, Reardon M, Gallagher MM, Hnatkova K, Camrn AJ, Malik M Circadiein Variation of the QT 
Interval in Patients With Sudden. Cardiac Death After Myocardial infarction Am J Cardiol 1 998, 81 .950-956, and Horns 
E, Marti V,- Guindb J, Laguna R Vinolas X, Caminal P, Elosua R, Bayes de Luna A. Automatic Measurement of corrected 
QT interval in Holter recordings; Comparison of its. dynamic behaviour in patients after myocardial infarction with and 
without life-threatening arrhythmias Am Heart J 1997;1 34.181 -7 However, theclinical value of many approaches to 
the assessment of the repolarisation dynamicity is unclear Reliable- descriptors to detect and quantify the autonomic 
modulations of ventricular repolarisation, would he extremely useful. The sensitivity. of VG and of the new decriptor 
TCRT for detecting autonomic modulations, could be translated into prognostic power for major cardiac events As 
mentioned above, the preliminary results indicate that such use of these descriptors is. certainly possible 
[0148] The conclusion from this study was that both VQ and TCRT are superior to QTmax and QTd in detecting 
autonomic modulations of myocardial repolarisation induced by postural changes and Nfelsalva manoeuvre The sen- 
sitivity of TCRT and the dynamicity of its reaction to. postural and autonomic modulations suggest that it would have 
particularly useful application for rateradaptive sensors in implantable antiarrhythmic devices Given its measurement 
could be simplified and made possible, from itracardiac leads it could be used for automatic detection of autonomic 
modulation of ventricular, repolarisation or other potentially arrythmogenic factors in implantable devices 
[0149] There now follows a description of preferred implantable devices, such as a pacemaker (see Figures 12 to 
14) and a pacemaker cardioverter-defibrillator (see Figures 15 and 16), which could use TCRT to improve pacing to 
act as a monitor for detecting ventricular repolarisation abnormalities or even act as an alarm to warn of autonomic 
conditions that pose a risk to the patient 

[0150] Figure 12 is a simplified schematic view of one embodiment of an implantable medical device ("IMD") 10 of 
the present invention IMD 10 shown in Figure 12 is a pacemaker comprising at least one of pacing and sensing leads 
16 and 18 attached to hermetically sealed enclosure 14 and implanted near human or mammalian heart 8 Pacing and 
sensing leads 16 and 18 sense electrical signals attendant to the depolarization and re-polarization of the heart 8, and 
further provide pacing pulses for causing depolarization ofcardiac tissue in the vicinity of the distal, ends, thereof Leads 
16 and 18 may have unipolar or bipolar electrodes disposed thereon, as is well known in the art Examples of IMD 10 
include Implantable cardiac pacemakers disclosed in US Patent No. 5,158,078 to Bennett et al , U.S. Patent No. 
5,312,453 to Shelton et al or U S Patent No. 5,144,949 to Olson, all hereby incorporated by reference herein, each 
in its respective entirety 

[0151] Figure 13 shows connector header module 12and hermetically scaled enclosure 14-of IMD 10 located in and 
near human or mammalian heart 8 Atria! and ventricular pacing leads 1 6 and 1 8 extend from connector header module 
1 2 to the right atrium and ventricle, respectively, of heart 8 Atrial electrodes 20 and 21 disposed at the distal end of 
atrial pacing lead 16 are located in the right atrium Ventricular electrodes 28 and 29 at the distal end of ventricular 
pacing lead 18 are located in the right ventricle 

[0152] Figure 14 shows a block diagram illustrating the constituent components of IMD 10 in accordance with one 
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embodiment of the present invention, where IMD 10 is pacemaker having a microprocessor-based architecture. IMD 
10 is shown as including activity sensor or. accelerometer 11. which is preferably a piezpcerarnicacceleromeler bonded 
to a hybrid circuit located inside enclosure 14. Activity sensor 11 typically (although not necessarily) provides a sensor 
output that varies as a function of a measured parameter relating to a patient's metabolic requirements Additional 
information based oh the patient's autonomic tone .could be provided here to improve the pacing. For the sake of 
convenience, IMD 10 in Figure 14 is shown with lead 18 only connected- thereto; similar circuitry.and.cpnnections.not 
explicitly shown in Figure-1.4 apply to lead 16. . 

[0153] IMD 10 in Figure 14is.most : preferably programmable by means of an external programming unit (not .shown 
in the Figures). One'such prograrnmer.is the commercially available Medtronic Model. 979Q- prbgrammeri which is 
microprocessor-based. and provides a series of encoded." signals to IMD 10,. typically thrbugh;a programming head 
which transmits or telemeters radio-frequency (RF) encoded signals. to. IMD. 1 0. Such a telemetry system is described 
inU S Patent No. 5, 31 2,453 to Wyborny etah, hereby incorporated by reference herein in its entirety The programming 
methodology disclosed in Wyborny et al.'s '453 patent is identified herein for illustrative purposes only Any of a number 
of suitable programming and telemetry methodologies known in the art- may.be employed so long as : the desired infor- 
mation is transmitted to and from the pacemaker. ' 

[0154] As shown in Figure 1 4, lead 18 is coupled to node 50 in IMD TO through input capacitor 52 Activity sensor 
-or accelerometer 1 1 is most preferably attached to a hybrid circuit located inside hermetically, sealed enclosure 1 4 of 
£| MD 1 6 The output signal provided by activity sensor 1 1 is coupled to jnput/output circuit 54 : 1 npuVoutput circuit 54 
^contains analog circuits for .interfacing to heart 8, activity sensor 11, antenna 56 and circuits for the,application of 
•stimulating pulses to heart 8. The rate of .heart 8 is controlled by sottware-implemented algorithms stored in microcom- 
puter circuit 58 

-[0155] . Microcomputer circuit 58 preferably comprises onboard circuit 60 and off-board circuit 62 Circuit 58 may 
^correspond- to a microcomputer circuit disclosed- in U.S Patent No. 5,312,453 to Shelton et al , hereby incorporated 
by reference herein- in. its entirety ' On-board circuit 60 preferably includes microprocessor. 64, system clock .circuit 66 
and oh -board RAM 68 and ROM 70 Off-board circuit 62 preferably comprises a RAM/ROM unit. On-board circuit 60 
and off-board circuit 62 are each coupled by data communication bus. 72 to digital . control lerAimer circuit 74 . Micro- 
computer circuit 58 may. comprise acustom integrated circuit device augmented by standard RAM/ROM components 
[01 56] Electrical, components shown in Figure .1 4 are powered by an appropriate implantable battery power source 
76 in accordance with common practice in the art For the sake of clarity, the coupling of battery powerto the various 
components of I MD 1 0 is not shown in the Figures Antenna 56 is connected to input/output circuit 54 to permit uplink/ 
downlink telemetry through RF transmitter and receiver, telemetry unit 78.. By: way of example, telemetry unit .78 may 
correspond to that disclosed in U S Patent No. 4,566,063 issued to Thompson et ai ..hereby incorporated by reference, 
herein in its entirety, or to that disclosed in the above-referenced '453 patent to Wyborny etal It is generally preferred 
that the particular programming and telemetry scheme selected permit the entry and storage of cardiac rate+esponse 
parameters. The specific embodiments of antenna 56, input/output circuit 54 and telemetry unit 78 presented herein 
are shown for illustrative purposes only, and are not intended to limit the scope of the present invention. 
[0157] Continuing to refer to Figure 14, V^p-and Bias circuit 82 most preferably generates, stable voltage reference 
and bias currents for analog circuits included in input/output circuit 54 Analog-to-digital converter (ADC) and multiplexer 
unit 84 digitizes analog signals and voltages to provide "real-time 1 telemetry intracardiac signals and : battery . endof- 
[ife (EOL) replacement functions. Operating, commands for controlling the timing of IMD 10 are coupled by data bus 
72 to digital controller/timer circuit 74, where digital timers and counters establish .the overall escape interval of the 
IMD 10 as well as various refractory, blahkirrg and other timing windows lor controlling the. operation of peripheral 
components disposed within input/output circuit 54. 

[0158] Digital controller/timer circuit 74 is preferably coupled to sensing circuitry, including sense amplifier 88, peak 
sense and threshold. measurement unit 90 and comparator/threshold detector 92 ..Circuit 74- is further preferably cou- 
pled to electrogram (EGM) amplifier 94 for receiving amplified and processed signals sensed by lead 18 Sense amplifier 
88 amplifies sensed electrical cardiac signals and provides an- amplified. signal to peak sense and threshold measure- 
ment circuitry 90, which in turn provides an indication of peak sensed voltages and.measured sense amplifier threshold, 
voltages on multiple conductor signal path 67 to digital controller/timer circuit 74. An amplified sense amplifier signal 
is then, provided to comparator/threshold detector 92. By way of example, sense amplifier 88 may correspond to that 
disclosed in U S Patent No. 4,379,459 to Stein, hereby incorporated by reference herein in its entirety. 
[0159] The electrogram signal provided by EGM amplifier 94 is employed when IMD 10 is being interrogated by an 
external programmer to transmit a .representatiqo of a cardiac analog electrogram. See, for example, U.S Patent No. 
4,556,063 to Thompson et al,, hereby incorporated by reference herein in its entirety. Output pulse generator 96 pro- 
vides pacing stimuli to patient's heart 8 through coupling capacitor 98 in response to a pacing trigger signal provided 
by digital controller/timer circuit 74 each time the escape interval times but, an externally transmitted pacing command 
is received or in response to other stored commands as is well known in the pacing art By way of example, output 
amplifier 96 may correspond generally to an output amplifier disclosed in U.S Patent No. 4,476,868 to Thompson, 
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hereby incorporated by reference herein in Its entirety. 

[01 60] The specific embodiments of input amplifier 88, output amplifier (output pulse generator) 96 and EGM amplifier 
94 identified herein are presented lor illustrative purposes only, and are not intended to be limiting in respect of the 
scope of the present invention. The specific embodiments of such circuits may not be critical to practicing some em- 
bodiments of the present invention so -Jong as they provide means for generating a stimulating pulse and are capable 
of providing signals indicative 

[0161] In some preferred embodiments of the present invention, IMD 10 may operate in.various non-rate-responsive 
' modes, including, but not. limited to, ; DDD, DDI , Wl , VOO. and WT: modes. In other, preferred embodiments of the 
present invention,- IMD 10 may operate in various .rate-responsive, including, but not limited to/DDDR, DDIR, WIR, 
- VOOR and WTR modes. Some embodiments of the, present invention are capable of operating in both .non-rate- 
responsive and. rate responsive modes. Moreover, in various embodiments of the present invention IMD10 may be 
programmably configured to operate sp that it varies the rate at which it delivers stimulating pulses to heart 8 only in 
response to one or more selected sensor outputs being generated One such response could be generated through 
measurement of TCRT This could be in response to changes in the autonomic tone of the patient or changes in direction 
of the repolarfsation waveiront Numerous pacemaker features and functions not explicitly mentioned herein may be 
incorporated into IMD 10 while remaining within the scope of the present invention 

[01 62] The present invention is not limited in scope to single-sensor or dual-sensor pacemakers, and is not limited 
to IMD's cqmprisrig activity prpressure sensors only. Nor is the present invention limited in. seope : to single-chamber 
pacemakers, single-chamber leads for pacemakers or single-sensor or dual-sensor leads for pacemakers. Thus, var- 
ious embodiments of the present invention may be practiced in conjunction with more than two leads or with multiple- 
chamber pacemakers, for example At least some embodiments of the present invention may be applied equally well 
in the contexts of single-, dual-, triple- or quadruple- chamber pacemakers or other types of. IMD's, See, for example, 
U.S. Patent No. 5,800,465 to Thompspo et al., hereby incorporated by reference herein; in. its entirety as are all U S 
Patents referenced therein 

[0163] IMD 10 may also be a pacemaker-cardioverter- defibrillator ("PCD") corresponding to any .of numerous com- 
mercially available implantable PCD's Various embodiments of the present invention may be practiced in conjunction 
with PCD's such as those disclosed in U/S Patent No, 5,545,186 to Ojson et al., U S Patent No 5,354,316 to Keim el, 
U SI. Patent No. 5,314,430 to Bardy, U S Patent No: 5,131,388 to P.less and U.S. Patent No. 4, .82.1 ,723 to Baker et 
ak, all hereby iricorporatied by reference herein, each in its respective entirety 

[01 64] Figures 1 5. and 1 6 illustrate one embodiment of IMD 10 and a corresponding lead set of the present invention, 
where IMD 10 isa.PCD. In Figure 15, the ventricular leadtakes the form of leads disclosed in U S. Patent Nos 5,099,838 
and 5,314,430. to Bardy, and includes an elongated insulatiyc lead body 1 carrying three concentric coiled conductors 
separated from, one another by tubular insulative sheaths Located adjacent the distal end pMead I are ring electrode 
2, extendable helix electrode 3 mounted retractably within insulative electrode head 4 and, elongated coil electrode 5 
Each of the electrodes is coupled to one of the coiled conductors within lead body t Electrodes 2 and 3 are; employed 
for cardiac pacing and for sensing ventricular depolarizations. At the proximal end of the lead is bifurcated connector 
6 which carries three electrical connectors, each coupled to one of the coiled conductors. Defibrillation electrode (elon- 
gated coil electrode) 5 may be fabricated from platinum, platinum alloy or other materials known to be usable in im- 
plantable defibrillation electrodes and may be about 1 5 cm in length. 

[0165] The atrial/SVC lead shown in Figure 15. includes elongated insulative lead body 7 carrying three concentric 
coiled conductors separated from one another by tubular insulative sheaths corresponding to the structure of the ven- 
tricular lead Located adjacent the J-shaped distal end of the lead are ring electrode 9 and extendable helix electrode 
1 3 mounted retractably within an insulative electrode head 15. Each of the electrodes is coupled to one of the coiled 
conductors within lead body 7 Electrodes 1 3 and 9 are employed for atrial pacing andf or sensing atrial depolarizations 
Elongatecf coil electrode 19 is provided proximal to electrode 9 and coupled to the third conductor within lead body 7 
Electrode 19 preferably is 10 cm in length or greater and is configured to extend from the SVC toward the tricuspid 
valve In one embodiment of the present invention, approximately 5 cm of the right atrium/SVG electrode, is focated in 
the right atnum with the remaining 5 cm located in the SVC At the proximal.end of the lead is bifurcated connector 17 
carrying three electrical connectors, each coupled to one of the coiled conductors. 

[0166] The coronary sinus lead shown in Figure 15 aissumes the lorm .of a coronary sinus lead disclosed in the above 
cited '838 patent issued to Bardy, and includes elongated insulative lead body 41 carrying one coiled conductor coupled 
to an elongated coiled defibrillation electrode 22 Electrode 22, illustrated in broken outline in Figure. 15, is located 
within the coronary sinus and great vein of theheart At the proximal end of the lead is connector. plug 23 carrying an 
electncal connector 41 coupled to the coiled conductor. The coronary sinus/great vein electrode 22 may be about 5 
cm in length 

[0167] Implantable PCD 10 is shown in Figure 15 in combination with leads 1, 7 and 41, and lead connector assem- 
blies 23, 17 and 6 inserted into connector block 12. Optionally, insulation of the outward facing portion of housing 14 
of PCD 1 0 may be provided using a plastic coating such as parylene orsilicone rubber, as is employed in some unipolar 
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cardiac pacemakers: The outward lacing portion,, however, may ^ be. left uninsulated-or some other division between 
insulated and uninsulated portions may be employed: The uninsulated portion of housing 14 serves as a subcutaneous 
defibrillation electrode to defibrillate either the atria or ventricles. Lead configurations other that those shown in Figure 
15 may be practiced in conjunction with the present invention, such as those shown in U S Patent No. 5,690.686 to 
5 Min etal.\ hereby incorporated. by reference herein in its entirety. . 

[0168] Figure 1 6 is a functional schematic diagram of one embodiment of implantable PCD 10 of thepresent invenr 
tion This diagram should .be takenas exemplary of the type of device in which various embodiments 61 the present 
invention may be embodied, and not as limiting, as.it is bejievedthat the.invention may be practiced in a wide vanety 
• of device implementations, including cardioverter and defibrillators which do not provide antMachycardia pacing ther- 

io. apies. . . 

roi69] IMD. 10 is provided with an electrode system: If the electrode configuration of Figure 15 is employed, the 
correspondence to the illustrated electrodes is as follows Electrode 25 in Figure .16 includes the uninsulated portion 
of the housing of PCD 1 0 Electrodes 25, 1 3/.21 and 5 are coupled to high voltage output circuit 27, which includes 
high voltage switches controlled by CV/defib control logic 30 via control bus 3.1 Switches disposed within circuit 27 

15 determine which electrodes are employed and which electrodes are coupled to the positive and negative terminals of 
the capacitor bank (which includes capacitors 33 and 35) during delivery of defibrillation pulses. 
[0170] Electrodes 2 and 3 are locatedon or in the ventricle and are coupled to the R-wave amplifier 37, which 
preferably, takes, the form of an automatic gain controlled amplifier providing an adjustable sensing threshold as a 
function of the measured R-wave amplitude. A signal is generated on R-out line 39 whenever the.signal sensed. bet ween 

20 electrodes 2 and 3 exceeds the present sensing threshold, • : 

[01 71] Electrodes 9 and 1 3 are located on or in the atrium and are coupled tp the P-wave amplifier 43, which preferably 
also takes the form of an automatic gain controlled amplifier providing an adjustable sensing threshold as a function 
of the measured P-wave amplitude . A signal, is generated on P^ut line 45 whenever the signal sensed between elec- 
trodes 9 and 13 exceeds the present sensing threshold. The general operation of.R-wave and P : wave amplifiers 37 

25 and 43 may correspond to that disclosed in U S. Pat. No. 5,11 7,824,. by Ke.mel et al. , issued Jun: 2, 1 992, for An 
Apparatus for Monitoring Electrical Physiologic Signals', hereby incorporated by reference herein in its ^rety 
[0172] Switch matrix 47 is used to select which of the available electrodes are coupled to wide band (0 5-200 Hz) 
amplifier 49 for use in digital signal analysis Selection ol electrodes .s controlled by the microprocessor 51 via data/ 
address bus 53, which selections may be varied as desired. Signals from the electrodes selected for coupling to band- 

30 pass amplifier (wide band amplifier) 49 are provided to multiplexer 55, and thereafter converted to multi-bit dig.tal 
signals by A/D converter 57, for storage in random access memory 59 under control of direct memory access circuit 
61 Microprocessor 51 may employ digital signal analysis techniques to characterize the digitized signals stored .n 
randomaccess memory 59 to recognize and classify the patient's heart rhythm employing any of the numerous signal 
processing methodologies known to the art For example. T.CRT could be. calculated and output signals generated 

35 rot73] m The remainder of the circuitry is dedicated to the provision of cardiac pacing, cardioversion and defibrillation 
therapies and for purposes^ the present invention may correspond to circuitry known to those skilled m the art The 
following exemplary apparatus is disclosed for accomplishing pacing, cardioversion and.defibrillation functions Pacer 
timing/control circuitry 63 preferably includes programmable digital counters which control the basic time intervals 

40 associated with DDD, Wl, DVI, VDD, AM, DDI and other modes of single and dual chamber pacing well known to the 
art. Circuitry 63 also preferably controls escape intervals.associated with anti-tachyarrhythmia pacing in bpththe atr.um 
and the ventricle, employing any anti-tachyarrhythmia pacing therapies known to the art. 

[0174] Intervals defined by pacing circuitiy 63 include atrial and ventricular pacing escape intervals, the refractory 
periods dunng which sensed P-waves and R-waves.are ineffective to restart timing of Ihe escape intervals and the 

45 pulse widths of the pacing pulses The durations of these intervals are determined, by microprocessor 51 , in response 
to stored data in memory 59 and are communicated to pacing circu itry 63 via address/data bus 53. Pacer circuitry 63 
also determines the amplitude o1 the cardiac pacing pulses under control of microprocessor 51 
[0175] Dunng pacing, escape interval punters within pacer timing/control circuitry 63 are reset upon sensing ol n- 
waves and P-waves as indicated by a signals on lines 39 and 45, and in accordance withlhe selected mode of pacing 

so on time-out trigger generation^ pacing pulses by pacer output circuitry 65a and 65b. which are. coupled to electrodes 
9 13 2 and 3 Escape interval counters are also reset on generation of pacing pulses and thereby control the basic 
timing of cardec pacing functions, including anti-tachyarrhythmia pacing The durations of the intewals defined by 
escape interval timers are determined by microprocessor 51 viadata/address bus 53, The value of .the count present, 
in the escape interval counters, when reset by sensed R-waves and P.waves may be used to. measure the durations 

55 of R-R inteivals, P-P intervals, P-R intervals, and R-P intervals, which measurements are stored in memory 59 and 
used to detect the presence of tachyarrhythmias. 

[0176] Microprocessor 51 most preferably operates as an interrupt driven device, and is responsive- to interrupts 
from pacertiming/control circuitry 63 corresponding to the occurrence sensed P-waves and R-waves and correspond- 
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ing to the generation of cardiac pacmg.pulses. Those interrupts are provided via data/address bus 53. Any necessary 
mathematical calculations to be performed by microprocessor 51 and any updating of the values or intervals controlled 
by pacer timing/control circuitry 63 take place following such interrupts 

[0177] Detection of atrial or ventricular tachyarrhythmias, as employed in the present invention, may correspond to 
5 tachyarrhythmia, detection. algorithms such as TCRT. Conventionally the presence, of an atrial . or ventricular tachyar- 
' rhythmia is confirrhed by detecting a sustained series of short R-R or P-P intervals of an average rate indicative of 
tachyarrhythmia or an unbroken series of short R^R or P-P intervals. The suddenness, of onset of the detected high 
urates, the stability of the high rates, and a number of other factors known in the artmay also be measured l.at this time. 
Appropriate ventricular tachyarrhythmia - detection. methodologies measuring, such factors are described, in U.S. Pat.. 
io : ,No. 4,726,380 issued to Vollmann, U:S. Pat. No. 4,880,005 issued to Pless et a! and U S Pat No. 4,830,006 issued 
to Haiuska et at.; all incorporated by reference . herein, each in its respective entirety. An additional set of tachycardia 
recognition methodologies is disclosed in the article "Onset and Stability for Ventricular Tachyarrhythmia Detection in 
ah Implantable Pacer-CardJOve^ter-Defibri.llator ,, by Olson etal. , published in Computers in Cardiology, Oct 7-10, 1986, 
IEEE Computer Society Press, pages 167-170, also incorporated by reference herein in its entirety Atrial fibnllation 
detection methodologies are disclosed in Published PCT Application Ser. No. US92/02829, Publication No. 
W092/18198, by Adams et at, and in the article "Automatic Tachycardia Recognition', by Arzbaecher eta!,, published 
in PACE, May-June, 1984, pp: 541-547, both of which are incorporated by reference herein in their entireties. 
[01 78] In the event an atrial or ventricular tachyarrhythmia is detected and an antkachyarmythmia. pacing regimen 
is desired, appropriate timing intervals for controlling generation of anti-tachyarrhythrnia pacing therapies .are loaded 
20 from microprocessor 51 into the pacer timing and control circuitry 63, to control the operation of the escape interval 
counters therein and to define refractory periods during which detection of R-waves and P-waves is ineffective to restart 
the escape interval counters 

[0179] • "Alternatively, .circuitry for. controlling the.timing and generation of anti-tachycardia pacing pulses as described 
in (J S. Pat No! 4,577,633, issued to Berkovits et al. on Mar. 25, 1986, U S. Pat! No. 4,880,005, issuedto Pless et al 
2S on Nov. 14 s 1989, U.S. Pat. No. 4,726,380 t issued to Vollmann ,et al. on Feb 23, 1 988 and U S, Pat. No. 4,5B7 ; 970, 
issued to Holley et al. on May 1 3, 1986, all of which are incorporated herein by reference in their entireties, may also 
be employed 

[0180] In the event that generation of a cardioversion or defibrillation pulse is required, microprocessor 51 may em- 
ploy an escape interval counter to control timing of such cardioversion and defibrillation pulses, as well as associated 

30 refractory periods. In response to the detection of atrial or ventricular fibrillation or tachyarrhythmia requinng a cardi- 
oversion, pulse, microprocessor 51 activates cardipversion/defibrillation control circuitry 30, which initiates charging of 
the high voltage capacitors 33. and 35 via charging circuit 69, under the control of high voltage charging control line 
71 : The voltage. on the high voltage capacitors is monitored via VCAP line 73, which is passed through multiplexer 55 
and in response to reaching a predetermined value set by microprocessor 51, results in generation of a logic signal 

35 on Gap Full (CF) line 77 to terminate charging Thereafter, timing of the delivery of the defibrillation or cardioversion 
pulse is controlled by pacer timing/control circuitry 63. Following delivery, ot the fibrillation or tachycardia therapy mi- 
croprocessor 51 returns the device to q cardiac pacing mode and awaits the next successive interrupt due to pacing 
or the occurrence of a sensed atrial or ventricular depolarization 

[01 8.1] S Weral. embodiments of appropriate systems for the.delivery and synchronization of ventricular cardioversion 

40 and defibrillation pulses and for controlling the timing functions, related, to them are disclosed in U S Patent No. 
5,188, 105 to Keimel, U S. Pat No, 5,269,298 to Adams et al and Li S Pat. No. 4,316,472 to Mrowski et al , hereby 
incorporated by reference herein, each in its respective entirety Any known cardioversion or defibrillation pulse control 
circuitry is believed to be usable in conjunction with various embodiments of the present invention, however. For ex- 
ample, circuitry controlling the timing and generation of cardioversion and defibrillation pulses such as that disclosed 

45 in U S Patent Nn 4,384,585. to Zipes, U S; Patent No. 4,949,719 to: Pless et at, or U S Patent No. 4,375,81 7 to Engle 
et al , all hereby incorporated by reference herein in their entireties, may also be employed 
[0182] Continuing to refer to Figure 1 6, delivery, of cardioversion or defibrillation pulses is. accomplished by output 
circuit 27 under the control of control, cir.cu.itry .30 via control bus 31 Output circuit 27 determines whether a monophasic 
or biphasic pulse is delivered, the polarity of the electrodes and which electrodes are involved in delivery of the pulse 

50 Output circuit 27 also includes high voltage switches which control whether electrodes are coupled together during 
delivery of the pulse Alternatively, electrodes intended to be coupled together during the pulse may simply be perma- 
nently coupled to one another, either exterior to or interior of the device housing, and polarity may similarly be pre-set, 
as in current implantable defibrillators. An example of output circuitry for delivery of biphasic pulse regimens to multiple 
electrode systems may be found in the above cited patent issued to Mehra and in U S. Patent No. 4,727,877, hereby 

ss incorporated by reference herein in its entirety, 

[01 83] An example of circuitry which may be used to control delivery of monophasic pulses is disclosed in U S.. Patent 
No. 5,163,427 to Keimel, also incorporated by reference herein in its entirety Output control circuitry similar to that 
disclosed in U S. Patent No. 4,953,551 to Mehra et al or U S Patent No. 4,800 ; 883 to Winstrom, both incorporated 
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by reference herein in their entireties, may also be used in conjunction with various embodiments of the present in- 
vention to deliver biphasic pulses. 

[0184] Alternatively, IMD 10 may be an implantable nerve stimulator or muscle stimulator such as that disclosed in 
U S. Patent No. 5,1 99,428 to Obel ef al., U.S Patent No. 5.207,218 to Carpentier et al or U.S. Patent No. 5,330,507 

5 . to.Schwartz, or an implantable monitoring device such as that disclosed in U.S. Patent No. 5,331 ,966 issued to Bennet 
: ' et"af.; all cif which ar^hereby incorporated by reference herein, each in-its respective entirety. The present invention 
is believed to find wide application toany form of implantable electrical device for.use in conjunction with electrical leads. 
[0185] As shown in Figures- 12 to 16, conventional implantable devices. 10 have several electrodes that could 'be 
. ■ . .■ used-for recording depolarisatibn and- repoiarisation signals for the purposes of deterniinirig TCRT or other descriptor 

10 • value. For. example; ih.the IMD.10 of Figure 1 3 signals could= be recorded between the atrial tip electrode 21 and the 
. ventricular tip electrode 29, between. the atrial tip electrode 2.1 and the IMD ("can") TO and between the ventricular tip 
electrode 29 and the IMD10 Other electrodes could be placed on the IMD 10, or within the heart 8 or remote from 
the heart connected via additional leads, The signals recorded in three or more directions could be processed to cal- 
culate vectors describing the. progress of the depolarisation and repoiarisation wavef rants. Microprocessor- 51 could 

is he used to transform the vectors into an optimum, dimensional space and descriptors such as TCRT and TMD could 
be determined. More complicated I MDs, for example, a pacemaker-cardioverter<lefibrillator, would have an additional 
lead or leads from which further data could be extracted. The a lnSync° implantable cardio defibrillator (ICD) which is 
sold by Medtfdriic, could record data between any two points of the ICD, SVC, three electrodes in the right ventricle, 
the leftventri&e electrode and two atrial electrodes. It is desirable to extract the wavefronl data from existing leads to 

20 avoid further obstructions beingintroduced into the heart pathways, for example, through the. presence of extra leads 
or thicker leads. In addition, the fewer the number of components, the less likely it is that the device will go wrong 
[01 85] In another embodiment, additional electrodes could be positioned on the patient corresponding, with- the or 
some:6f the standard independent ECG leads Information from these be fed to the IMD 10 for controlling pacing 
in response to autonomic changes or the detection of ventricular repoiarisation. abnormalities, in a f urther embodiment 

25 the leads of an IMD 1 0 could be used to provide data which is fed to external monitoring equipment. This might be in 
the form of a possible device which the patient wears so that the patient can be warned when TCRT values are at 
dangerous levels as a. result of autonomic changes, for example,, during exercise 

[0187] In another embodiment, the IMD 10 may record the data with respect to time which is then fed to an external 
data processing system for analysis, for example, during a check up. Such data may provide eairly indications of heart 
30 problems in addition to that being treated by the implantable device 10 These .might include ischemia and other con- 
ditions where damage to the heart muscle is caused which affect the direction of the repoiarisation wavefront, change- 
the autonomic tone orthe response ofttie autonomic system. : 

[0188] The descriptor TCRT could be calculated and used by itself in such situations or the ECG data could be 
processed further to extract additional descriptor values, which together withTCRT may. provide a better diagnosis tool 

35 [0189] In one embodiment envisaged in the p resent invention, a descriptor such as TCRT or TMD, or a combination, 
of descriptors, may he used to/control a drug delivery mechanism This may he by way of an implantable device or 
possibly via external equipment, for example, when the patient is in an intensive c^re unit. Such drugs may be used 
to reduce, clotting of blood (e.g asprin) or perhaps reduce the stress.levels in the patient. Sympathetic and parasym- 
pathetic agents could also be administered Where no correction of symptoms is observed, alarms could be triggered 

40 and assistance summoned . - 

[0190] As mentioned above,. TCRT has been found to provide a. useful measure of the autonomic tone in the patient. 
Nearly all anaesthetic agents create a sympathetic or a. parasympathetic effect in the autonomic system of a patient 
It would therefore be useful for the autonomic response, of the patient to be tested prior, to being anethetised. TCRT 
could be monitored while a patient performs Valsalva manoeuvre (or some other standard test) to check the autonomic 

45 response, prior to being anethetised for surgery. 

[0191] A third case study was carried to investigate whether QT dispersion could represent properly interlead het- 
erogeneity of ventricular repoiarisation The results obtained were, compared to. determine, the extent that the non- 
dipolar components of an ECG differ from the dipolar components, i e the residual energy of the T wave was investi- 
gated. The study was as follows. 

so [01 92] The concept ofthe so-cal led QT dispersion has recently attracted significant attention from the.clinical research 
community Various methods have been proposed to evaluate QT dispersion from the standard 12 lead electrocardi- 
ogram (EGG) by measuring QT intervals in individual leads Most frequently, the simple range of the QT interval meas- 
urements is used There -are numerous studies- indicating the-clinicai value of QT . dispersion Among others, increased 
QT dispersion has been reported to be associated with QT interval prolongation due to drugs of known proarrhythmic 

55 properties and to be less increased on drugs with lesser proarrhythmic effects, to predict mortality in general, epidemi- 
ological studies to identify patients who are at greater risk after surviving acute myocardial infarction and to mark 
therapeutic efficacy in the idiopathic long QT syndrome Recently, however, reports have also appeared challenging 
the clinical usefulness of QT dispersion 
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[0193] Since the introduction of the concept of QT dispersion, it has been speculated that the increased: range of- QT 
interval measurements is caused by the regional heterogeneity of the duration of ventricular repolarisation. It has been 
proposed that different leads. of the standard 12-lead ECG project the repolarisation signals of different regions of the 
myocardial tissue and that, consequently, increased dispersion is a sign of regional differences in the duration of re- 
5 polarisation. Indeed, studies comparing the QT dispersion with the. dispersion of these duration of monophasic action 
potentials found a geheraf correlation supporting this ■hypqth^fe/lrwas , bbs'erved;that QT dispersion is increasingly 
: prolonged with increasing.differences in the duration of monophasic action potentials, recorded at different endocardial 

' : sites. ■ • ■ • =■■■■•■ ' \\ : ■ : \- , . 

[01 94] At tM sametime, these studies do.not offer direct.proof.1hat. increased dispersion of the QT intervals in the 
io standard 12-lead ECGs : measures directly.the same phenomena as the dispersionof durations of monophasic. action, 
potentials. If an increased heterogeneity exists in thedurations of monophasic action potentials, the repolarisation 
sequence is more disturbed, the vectorcardjographic loop, of the Twave is more abnormal and the projections of this 
loop into the standard ECG leads are more complicated than in normal electrocardiograms. 
[01 95] Recent studieshave shown that a similar value of QT dispersion is recorded in full 1 2-fead ECGs and in their 

is reconstruction from orthogonal XYZ leads and that QT dispersion values correlate with parameters of the vectorcar- 
diographs T loop morphology. It has consequently been speculated that the different projections of the T wave vector 
onto the different leads of the standard ECG play an essential role, and that the hypothesis of QT dispersion repre- 
senting a direct measure of the heterogeneity of ventricular repolarisation durations is. flawed. Such a concept can 
explain even the earlier observations of the correlation of QT dispersion with the heterogeneity of monophasic action 

20 potential durations However, the studies of orthogonally reconstructed 12-lead ECGs and of correlations with T loop 
morphology only prove that the projections of the T wave vector play an .important role in determining QT dispersion 
but they do not prove that a regional heterogeneity of myocardial repolarisation. duration is not involved at all. , 
[0196] To solve the problem of whether QT dispersion is, or is not, associated with regional heterogeneity of myo- 
cardial repolarisation, potentially in addition to the T wave vector . projection, a direct study was conducted comparing 

25 qt dispersion with: electrocardiographic signals that are not attributable to the orthogonal vector of the T wave. In 
12-lead ECGs obtained from several clinically well defined populations, QT dispersion and, using special signal 
processing techniques, the extent of non-dipolar components: by- which the individuals leads of a 12-lead ECG differ 
from the 3^D vector of the repolarisation signals were measured. 
J0197] The study involved four separate groups of subjects. 

30 [0198] The group of normal subjects consisted of 78 normal healthy volunteers (aged 47 ± 16 years, 23 women) 
with normal physical examination and normal 12-lead ECG. At the time of the study, none of the normal subjects was 
on any medication and on the day of the study, the subjects were asked to refrain from smoking and from alcohol and 
caffeine intake 

[01 99} The group of hypertrophic cardiomyopathy (HCM) patients consisted of 68 patients (mean age 38 ± 1 5 years, 

35 21 women) referred to St George's Hospital London,. England, for diagnosis, risk stratification, management of symp- 
toms, and/or follow-up evaluation Following the established guidelines, the diagnosis of HCM was based on the pres- 
ence of left ventricular hypertrophy on 2D echocardiography in the absence of other cardiac or systemic disease that 
may cause left ventricular hypertrophy For ethical reasons; patients were not required to discontinue therapy before 
this study At the time of the study, 9 patients were on medication with established or potential effects on myocardial 

40 repolarisation (amiodarone, n = 6; and sotaloi, n = 3) 

[0200] The group of idiopathic dilated cardiomyopathy (DCM) patients consisted, of 72 patients. (mean age 48 ± 1 5 
years, 29 women). The diagnosis of idiopathic DCM was based on enlarged left. ventricular diameters (left ventricular 
diameters, diastolic 64 ± 10 mm, systolic 51 ± 13 mm) with reduced systolic function without any underlying causes 
of DCM At the time of the recording, 16 patients were on amiodarone 

45 [0201] Finally, the group of. patients with acute myocardial infarction (AMI) consisted of 81 patients (mean age 63 ± 
12 years, 20 women) Diagnosis of acute myocardial infarction was based on. previously. pubJished criteria, i.e. the 
presence of at least 2 of 3 standard signsof (a) typiqal chest pain, (b)cardiac enzymes elevated more than twice above 
to normal levels of pur laboratory and (c) typical ECG changes. History ol previous myocardial infarction was recorded 
in 1 3 patients, 44 patients had an anterior infarction, and at the time of hospital admission, 68 patients received thronv 

50 bolytic therapy At the time of the study,, none of the patients was on an antiarrhythmic therapy, 80 patients were receiving 
aspirin, 43 beta-blockers, 32 diuretics, and 34 ACE inhibitors. 

[0202] In the patient groups, subjects were not eligible for this study if in atnalfibrillation or other non-sinus, rhythm, 
in the presence of atrioventricular conduction block, or-with a QRS duration > 120 ms, 

[0203] in each subject of each group, 10 serial 12-lead ECGs were recorded in a supine resting position using a 
55 digital 12-lead electrocardiograph MAC VU by Marquette Medical Systems (Milwaukee, Wisconsin) Each ECG re- 
corded simultaneously ail 12-leads for 10 seconds and the serial ECGs were performed one after another without 
removing the electrodes In all subjects, all 10 ECGs were recorded within less than 3 min 

[0204] All ECGs were recorded after careful skin preparation The healthy subjects were recorded after being m- 



28 



structed to refrain from smoking and caffeine intake on the day of the study; HCM and.DCM patients were recorded at 
the time of presentation at a specialised out-patient clinic of our Hospital. The AMI patients were recorded on day 1 
following the index infarction 

[0205] Each ECG was stored on a floppy disc (500 Hz sampling at 1 2 bit resolution) and transferred to a dedicated 

5 . workstation equipped with the QT Guard package (Marquette Medical Systems) which was used to construct the so- 
•cair^' mediaabS'af of each lead of each elebtrocadidgram:-These : median-beats represent an ideal- QRST- complex 
of each lead of the ECG and, compared to the native ECG signal, have an improved signal to noise ratio.. The median 
beats ; were further used to measure the.QT dispersion and.thenon<lipolar. components ot each.ECG.- . .; 
. . [0206] In each electrocardiogram, QT dispersion was. measured using. the QT Guard package.- In each leadofeach 

ro -■' ECG the noise of the isoelectric, line was. measured and compared, with the voltage of the peak of the T wave. If the 
standard deviation of the T-P segment signal did not exceed 70% of the maximum T wave amplitude and if the T wave 
amplitude was > 60mV, the leadwas. measured, otherwise it was excluded from: the measurement !n each lead, the. 
end of T wave-was firstly determined automatically using the intersection of the isoelectric line with the tangent to the 
inflection point of the/descending part of. the T wave. (The tangent was calculated using least square fit-to the 3 samples 

is above and 3 samples below the inflection point. ) These automatic measurements were subjected to visual checks by 
an experienced and electrocardiographically trained operator who corrected the automatic measurements manually 
when necessary During the editing process, the operator was kept blinded with respect to the association of the ECGs 
L-.with the individual groups in -the study though clinical ECG diagnosis was. possible in a number of cases. 
40207] An ECG was accepted for the QT dispersion measurement if the T wave offset was. measured in at least 9 

20 l 0 f the 12 standard leads. In such a case,. QT dispersion was expressed using three different methods': as the range of 
the QT interval durations in alt measured leads (that is the difference between.the maximum andminimum QT interval 
measured - QTd method 1 ), as the standard deviation of the QT interval durations in all measurable leads (QTd method 
2) and as the difference between the upper and lower quartile of the QT interval durations in all measurable leads 
" (QTd method 3) Methods 2 and 3 for expressing QT dispersion were used in an attempt to overcome the technical 

ss problems associated with the simple measurement of QT interval range In addition to these measures of QT dispersion, 
maximum QT interval; was taken as the maximum of the QT intervals of all measurable leads. Heart rate was also 
derived from each ECG 

[0208] For each method of QT dispersion measurement, the . results obtained in the: serial electrocardiograms of the 
same subject wereaveraged and the mean value was used as the true measure of QT dispersion for the given individual 

30 The representative values of maximumQT interval and of heart rate were obtained for each subject in the same way. 
[0209] The concept attributing QT dispersion to the regional differences of myocardial repolansation assumes that 
in addition to the global T wave vector, each lead of the 1 2-l.ead ECG records signals from a region of the heart (nearest 
to the electrode in case of the precordial leads) which are not recorded by any other lead. To quantify the presence of 
such signals, we have measured the non-dipolar components of Ihe 12-lead ECG, i.e the extent of the residuum of 

35 the T wave which exists after subtracting the T wave vector. 

[0210] More specifically, using the technique described above, in relation to the first study, the signals from the eight 
independent leads of the 12-lead ECG (namely leads I, II, VI, V2, ,V6) were subjected to the Singular Value Decom- 
position and the electrocardiogram reconstructed in an orthogonal 8-lead system In such a system, the first lead 
contained the maximum energy in one single direction, the second lead the maximum energy perpendicular to the first 

40 lead the third lead, the maximum energy perpendicular Id the first two leads, etc In this way, the energy embedded 
in the first three orthogonal leads corresponded to the energy of the T wave vector while the energy in the remaining 
leads 4 - 8 corresponded to the non-dipolar components summed over all 1 2-leads of the onginal ECG (see Figure 2a 
which illustrates the recorded. ECG signals and Figure 2b which illustrates those signals when they are reconstructed 
using Singular Value Decomposition) The Singular Value Decomposition is dependent on the region of the ECG signal 

45 lor whibh the optimisation of the orthogonal leads is performed For the purposes of this study, we have optimised the 
Singular Value Decomposition for the T wave rather than for the QFIS complex 

[0211] For each ECG, the proportion between the non-dipolar components in orthogonal leads 4 - 8 and (that is the 
time integral of leads 4 - 8 within the T wave) and the energy of the T wave vector in leads 1 - .3 (that is thefime mtegral 
of leads 1 - 3 of the optimised orthogonal system of over the same time) was obtained Similar to the measurement of 
so QT dispersion, these values obtained from the serial ECGs in. each subject were averaged and the result taken as the 
true measure of the non-dipolar component. 

' [0212] For the puiposes of this study, we term the proportion between the non-dipolar and 3D vector components 
"th e relative T wave residu um". 

[0213] Subjects were excluded if at least 5 of the 10 serial ECGs either did not provide QT dispersion measurement 
ss based on the acceptance criteria as above, or were rejected by the singular value decomposition package because of 
low signal to noise ratio or other technical reasons. 

[021 4] The values of heart rate, Fridericia corrected, maximum QT interval (QTc). QT dispersion and. of the relative 
T wave residuum were compared in individual groups of the study. Since the distribution of the values of the relative 
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T wave residua is not known, thenon^aramet^ for this purpose. 

The correspondence between QT dispersion and T wave residua was ; examined using Spearman rank correlation 
coefficients which were calculated for the complete study as well as for the individual clinically defined populations In 
the same way, the correspondence between T wave residua.and heart rate and QTc interval were evaluated. 
[0215] Unless specified otherwise,, the data-in tables are presented as mean ± standard deviation while in whisker 
charts, data are presented as mean ± standard error of the mean. A p value less than 0 05 was considered statistically 
significant. *■'•*'. •./ ' . : * • 

[021.6] For ECG processing reasons (mainly for f waves of too low amplitude in too many leads), 5 HCM patients !. . 
• and 10.DCM patients were excluded from the analysis, in the remaining subjects (78 normal volunteers, 63 HCM . 
patients; 62 DCM patients, and 81 AMI patients), fewer than all serial ECGs were used in 3 DCM, 1 HCM, and 2 AMI 
patients 

[021 7] Table 7 shows the correlation coefficients between the individual methods for expressing QT dispersion While 
Method 1 (range) was very closely correlated' with Method 2 (standard deviation), Method 3 (inter-quartile difference) 
leads to a somewhat less close correlation although the relationship remains very strong and very statistically significant 
(Figures 17a and 17b) 

[0218] Table 8 shows the differences between the study populations in respect of heart rate and QTc. interval. Figures 
18a and 18b shows the values of QT dispersion (Method 1 and Method 3) in the four populations of the study. With 
Method 3 (inter-quartile- range), all the differences between individual pairs of populations were. statistically significant 
. or nearly statistically significant. However, Method 1 (range of QT intervals) and Method 2 (standard deviation) did not 
differentiate between normal subjects and DCM patients (p = 0,92 and p = 0 35 for Method 1 and Method 2, respec- 
tively) Note that the results for Method 2 of QT dispersion are not shown in the figure -. although numerically different, 
the: values practically reproduced the comparisonswith Method 1 . 
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Correlation, coefficients between QT dispersion indices 


Group 


QTd 1 vsQTd2 


QTd 1 vs QTd 3 


QTd 2 vs QTd 3 


NRM 


: 0 9824 


0 41.36 


0 4889 


HCM 


0 9689 


0 6677 


0 7725 


DGM 


0 9576 


6 6437 


0 7529 


AMI 


0.9655 


0.7772 


6.8685 


. Total population 


0 9812 


07212 


0 7942 



QTd. 1 = range of measurable. QT intervals,. QTd 2 standard deviation of measurable QT intervals, QTd 3 = inter- 
quartile difference of measurable QT intervals NRM = normal healthy volunteers, HCM = hypertrophic cardiomyopathy 
patients, DCM = idiopathic dilated cardiomyopathy patients, AMI = survivors of acute myocardial infarction 
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Comparisons between heart rate and QTc interval durations in study groups 


Group 


Heart rate (bpm) 


QTc interval (ms) 


NRM 


66.7 ±9 4 


406.5 ±1.7 7 


HCM 


67 6±12 9 


447 2 ±26 5 


DCM 


76 0 ±13.0 


429 4 ± 35.3 


AMI 


74 1+14.7 


441 1 ±33 8 



[021 9] All differences between heart rate in individual groups were statistically significant with the exception of normal 
subjects vs HCM patients, DCM vs AMI patients All differences between QTc interval were statistically significant with 
the exception of HCM vs AMI patients NRM = normal healthy volunteers, HCM = hypertrophic cardiomyopathy patients, 
DCM = idiopathic dilated cardiomyopathy patients, AMI = survivors of acute myocardial infarction bpm = beats per 
minute 

[0220] Figure 19 shows the values of the relative T wave residua in the individual populations of the study. With the 
exception of the difference between normal subjects and HCM patients (p = 0 1.4), ail the differences between individual 
populations were statistically significant. It should be noted that the values of the relative T wave residuum are very 
small, the mean value in normal healthy subjects being approximately 0 03% which means, that in the normal subjects, 
we have found the proportion between non -dipolar. and dipolar components of the T wave in the order of 3 in 1 0,000 
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[0221] The correlations between the relative T wave residua and the measures of QT dispersion is presented in 
Table 9. Corresponding scatter diagrams arc shown in Figures 20a and 20b In individual populations of the study, 
statistical significance of the correlation between relative T wave residua and.QT dispersion was only reached in HCM 
patients. Note that in the DCM and HCM populations, an opposite relationship between QT dispersion and relative T 
5. . wave residua was observed. Statistical, significance^ the correlation was reached in the total population of the study 
almost, certainly because ofsimilar trends from ^ ..■».✓■.■....■■ -.- 

[02.22] Correlation coefficients of the relative T wave residua with heart rate and maximum QTc interval are also 
-shown in Table 9. Similar to the measures of QT dispersion; the residuum is related neither .to heart rate nor to the QTc 
. interval ' 

10 [(5223] The findings of this third case study may be summarised as follows: 

a) The non-dipolar components, (i.e electrocardiographic regional heterogeneity) of the repolansation signals are 
• measurable in digital 1.2^lead ECGs. 

b) These nonrdipolar components. differ in different clinically well-defined groups. 

is cj The so-called QT dispersion is unrelated to the non-dipolar components of the T wave Consequently, QT dis- 

persion does not represent a direct measure Of regional heterogeneity of ventricular repolansation 
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Correlation coefficients between QT dispersion and relative T wave res 


duum 


Group 


QTd.1 vs Twr 


QTd2 vsTwr 


QTd3 vs Twr 


HR vs Twr 


QTc vs Twr 


NRM 


-0 0446 
NS 


-0 09.45 
NS 


-0 0811 
• NS 


0.0794 
NS 


P2193 . 
NS 


J I CM 


0 2805 
P = 0026 


0.2882. 
P = Q 022 


0 3305 
P = 0 0p8 


-0 2027 
NS 


0 1322 " 
NS 


DCM 


-0 1531 
NS 


-Q 1755 
NS 


-0 2201 
• NS 


-0 0873 
NS 


-0.0317 
NS 


AMI 


0 0771 
NS 


0 0445 
NS 


0 0393 
NS 


0 1054 

NS . 


0.2807 
P = 00.11 


total population 


0 2165 
p = OO0O26 


0 2380 
p ■= 6x1 0^ 5 


0:2982 
P= 3x10-7 


0.0135 
NS 


0.3270 
P = 2x10" 8 



QTd 1 = range of measurable QT intervals, QTd 2 = standard deviation of measurable QT intervals, QTd 3 = inter- 
quartile difference of measurable QWT intervals HP = heart rate, QTc = Fridericia corrected maximum QT interval 
Twr = relative T wave residuum. NRM = normal healthy volunteers, HCM = hypertrophic cardiomyopathy patients, 
DCM = idiopathic dilated cardiomyopathy patients, AMI = survivors of acute myocardial infarction^ 
[0224] The findings of this study shows that QT dispersion is largely caused by the different projection of the T wave 
vector in different leads of the standard ECG These findings aiso show a significant difference between QT dispersion 
in groups, with a different morphology of the vectorcardiograph^ T wave loop Also, this is the firsf study to show that 
in addition ib the projection effects 6i the J wave vector, no regional components of signal repolansation play a role in 
determining QT dispersion 

[0225] This case study does not disprove the clinical utility of QT dispersion Indeed, the large number of clinical 
studies showing the potential of QT dispersion are consistent with the hypothesis that patients at greater risk (e g 
patients on proarrhythmic therapy or patients with advanced ischaemic heart disease liable to ventricular tachycardia/ 
fibrillation) have a more complex T wave vector and theref ore a more complex projection of the T wave vector into the 
individual ECG leads Perhaps, practical aspects of QT dispersion measurement also play a role More complex pat- 
terns of the ECG repolarisation signals may clearly lead to increased difficulty with determining the end of the T wave 
and may consequently result in the measurement of an increased QT dispersion If this is the case, the more complex 
projections of the T wave vector arc combined with a systematic bias towards increased QT dispersion values in 
patients with disturbed ventricular repolarisation. The studies -reporting poor Intra- and r inter<)bserver reproducibility 
and poor intrasubject stability of QT dispersion, assessment point in this direction One study found not only increased 
QT dispersion in HCfyl patients compared to normals but also a lower reproducibility of QT dispersion measurement 
in these patients 

[0226] The study, together with previously published shows that QT dispersion may relate to nothing, more than an 
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expression of T wave loop abnormalities, rather than T wave loop morphology: Although technically more difficult to 
quantify^ the morphology of the T wave loop appears to be a far valuable ECG factor. 

[0227] In the study it was observed that the relative T wave residuum differs between different clinically well defined 
groups of our study. T wave. residuum may be useful clinically as a diagnostic tool, either on its own, or more preferably, 
in conjunction with other descriptors which could indicate certain heart conditions. It is possible to take the view that 
the relative t wave residuum truly ' corresponds to the local heterogeneity in ventricular repo^ 



/ Claims . . . r. 

1. A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals hedrt from different spatial positions with respect to the heart during depoiarisation and repolarisation of 
the patient's heart, the plurality of electrical signals monitoring the propagation of depoiarisation and repolarisation 
waves originating in. the patient's heart, processingthe plurality of electrical signals to yield a vector which describes 
the propagation direction of one of the depoiarisation and repolarisation wavefronts, and a set of a plurality of 
vectors which describe the propagation direction of the other of the, depoiarisation and repolarisation wavefronts 
with respect to time, and determining the vector deviation between the depoiarisation and repolarisation wavefronts 
by measuring the angle between pairs of. respective vectors for all combinations of depoiarisation vector to repo- 
larisation vector between predetermined time limits. 

2. A method as claimed in claim 1 , wherein the cosine of the angle between each depoiarisation / repolarisation 
vector pair is calculated. 

3. A method as claimed in claim i or 2, wherein angles are measured between a vector for the repolansation wave, 
which corresponds to a direction of maximum energy of the repolarisation wave, and each vector from a set of 
vectors describing, the depolarisatioh wave at different time .instances, the set of vectors corresponding to sub- 
stantially the whole duration of depoiarisation of the patient's heart 

4.. A method as claimed in claim 1, 2 or 3, wherein the mean cosine of the angle is calculated for all depoiarisation / 
repolarisation vector pairs. 

5. A method of characterising ventricular operation of a patient's heart, comprising sensing the propagation of depo- 
lansatipn.and repolarisation. waves originating in the heart, determining vectors which arc representative of the 
direction of the wavefronts of the depoiarisation and repolarisation waves, and determining the vector deviation 
between the depolarisatioh and repolarisation vectors by determining the cosine of the angle between the vectors 
describing the depoiarisation and repolarisaition wavefronts 

6. A method as claimed in claim 5, wherein the vector (deviation is a function of. 

a) the cosine of the angle between two vectors, each vector describing one of the depoiarisation and repolari- 
sation wavefronts, 

b) the cosines of the angles between a vector describing either the depolarisatioh or repolarisation wavef ront 
and a set of vectors describing the other of the depoiarisation or repolarisation wavefront for a plurality of time 
instances; or 

c) the cosines of the angles between a set of vectors describing the depoiarisation wavefront for a plurality of 
time instances and a set of vectors describing the repolarisation wavefront for a plurality of time instances. 

7. A method as claimed in claim. 5 or 6, wherein data is produced for the propagation of. the depoiarisation and 
repolarisation wave with respect to a first set of axes, the data is transformed to a new set of axes defining an 
optimised orthogonaldomain having a first axis aligned with a direction of maximum energy for depoiarisation or 
repolarisation and wherein the vector deviation is measured in the optimised orthogonal domain 

8. A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals from different spatial positions with respect to the heart during depoiarisation and repolarisation of the 
heart, the plurality of electrical signals being sensed by an implantable medical device, and being associated with 
the propagation of depoiarisation and repolarisation waves originating from a patient's heart, processing the elec- 
trical signals, to yield a plurality of vectors which describe the propagation directbn of a wavefront for a depoiari- 
sation wave and a plurality of vectors which describe the propagation direction of a wavefront for a repolarisation 
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wave, wherein ventricular operation is characterised. in terms of the cosine of the angle between the. plurality of 
vectors for the depolarisation and repoiarisation waves'. ' ■'•= ."^ •"; 

9 A method as claimed n cla.m 8. wherein at least one of the plurality of vectors describes the propagation of the 
' wavefront asa function of time and the mean of .the cosine of the angle between pairs of vectors is determined. 

10 An .implantable medical devbe comprising a plurality of medical, electrical leads, the leads having electrodes for 
■ V sensing electrical signals from different spatial positions in, on. or near a patient's heart,- wherein the device proc- . 
. esses the eiectrical signals to yield directions: of propagation for depolarisation and repolarisationwaves of apa- 

: : tient's heart, calculates the angle of deviation between the depolarisation and repolarisat.cn waves, and generates 
an output signal corresponding to the angle of deviation 

11. A device as claimed in cla.m 10. wherein the outputsignal varies ..n accordance with the cosine of the. angle of 
deviation 

12 A device as cla.med in claim 11 . wherein the angle of deviation between the depdlansation and ^polarisation 
. .. ." waves is calculated with respect to time and the output signal varies in. accordance with the mean, of the cosine 

of the angle of deviation for a plurality o.f time instances. 

13 A device as claimed in claim 1 2, wherein the output signal varies in accordance with the mean of a set of cosine 
■ ' valuesforthe angle of deviation between indirection of the repoiarisation wavefront for a maximum energy value 

and the direction of the depolarisation wavefront as afunctfen of timefor all time, instances between startand f.n.sh 
points of depolarisation. 

14. A device as claimed in any of claims 10to 13, wherein the output- signal is usedforat least.pne.of. 

a) control pacing of the patient's heart, 

b) monitor the condition of the patient's heart, 

cj monitor the progression of disease in the patient's heart;. 
30 d) raise an alarm, when the angle is outside predetermined limits, 

e.) control a drug dispensing, pump/and / or; 
f ) monitor the response of the patient's autonomic system, 

15 A device as claimed in any of claims 11 to 14, wherein the device is selected 1rom the group, of, pacemaker, cardi- 
35 ' overter, defibrillator, pacemaker-cardioverter-defibrillator, heart monitor and drug dispensing pump 

.16. A device as claimed in any of claims 11 to 15. wherem. additional electrical signals are sensed by: 

a) at least one electrode provided on a housing of the device, 
40 b) at least one subcutaneous electrode, and / or 

c) at least one external electrode that is applied to a patient's body. 

17. A device as claimed in claim 16, wherein the at least one electrode is coupled to the device by a medical electrical 
lead, by an electncal connection or by a radio frequency transmitter 
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18. A method of determining whether a patient is suffering from heart failure, and preferably ^Mhe^patient * 
suffering from one- of hypertropic cardiomyopathy, idiopathic dilated cardiomyopathy and acute myocardial infarc- 
tion, comprising .mea Jriruj the cosine of the angle between propagation directions of depolarisation and repoiari- 
sation waves produced by the patient's heart 

19 A method of determining whether a patient .s suffering from heart failure, and preferably whether fhe^ patient 
suffering from one of hypertropic cardiomyopathy, idiopathic dilated cardiomyopathy endnote 
tion comprising sensing a plurality of electrical. signals from. different, spatial posrt,ons in, on. .or near* patient s 
Tart aviated wKh. depolarisation andrepolarisation o, the patient's heart by emptying an imp 
device which is capable of sensing depolarisation and repoiarisation waves and processing the plurality -o signals 
to determine the angle between propagation directions of depolarisation and repolaneaton waves produced by 
the patient's heart 



33 



20. A method of monitoring changes in the autonomic system of a patient comprising sensing a plurality of electrical 
signals associated with the propagation direction of depolarisation and repolarization waves originating in the 
patient's heart, determining the cosine of the angle of deviation between the depolarisation and repolarisation 
waves, and monitoring changes in the cosine of the angle of deviation, wherein preferably said method is accom- 

• • . plished by employing an implantable medical device which is capable of. sensing the waves and processing.same. 

21. A method of characterising ventricular operation, comprising sensing a plurality of electrical signals associated 
with the propagation of a repolarisation wave originating in the. patient's heart, the. plurality of electrical signals , 

. being sensed from different spatial petitions on, in or near the. patient's heart, processing the plurality of electrical 
• signals to yield: a plurality of vectors that. are. representative of the wavefront of the repolarisation wave, and de- 
termining a measure of the spatial variation of the repolarisation wavefront. 

22. A method as claimed in claim 21, wherein the spatial variation is calculated by determining vector contributions 
for the repolarisation wavefront in. each of a set of predetermined directions and measuring the angle between 
pairs of vector contributions 

23. A method as claimed in clairri 22, wherein the pre-determined directions, correspond, with at least three of the 
standard ECG channels of I, II, VI , V2, V3, V4, V5 and V6, and preferably the vector contribution for the ECG 
channel of V1 Is ignored in the calculation 

24. A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals to monitor repolarisation of the heart from different spatial positions with respect to the patient's heart, 
processing the plurality of signals to yield a vector describing the propagation of a repolarisation wave, through the 
heart, projecting the vector onto. a -set of axes to determine vector contributions of the signal vector in the directions 
of the axes, and measuring the angle between pairs of vector contributions. 

25. A method, as claimed in claim 24, wherein the vector corresponds to a direction of maximum energy of the repo- 
larisation wave 

26. A method of characterising ventncular operation of a patient's heart comprising sensing a plurality of electrical 
signals to monitor propagation of repolarisation through the heart from different spatial positions with respect to 
the patient's heart, processing the plurality of electrical signals to yield a vector describing the propagation of a 
repolarisation wave with respect to time and with respect to a first set of axes defining an optimum domain space, 
mapping the path of a tip of the vector in the optimum domain space to generate a T-wave loop and calculating a 
parameter describing the morphology variation of that loop. 

27. A method as claimed in claim 26, wherein the parameter is determined by projecting the T-wave loop on to recon- 
struction vectors corresponding to electrode positions to generate vector contributions in those electrode directions, 
and determining the angle between all pairs of vector contributions. 

28. A method as claimed in claim 27, wherein the reconstruction vectors correspond to the position of at least three 
of the standard ECG leads of I, If, V1, V2, V3, V4, V5 and V6, and preferably the vector contribution, from the 
standard ECG lead of V1 is. ignored in the calculation 

29. A method as claimed in any of claims 26 to 28, wherein the T-wave loop is mapped in two orthogonal dimensions 
and the energy of the T-wave loop in the two orthogonal dimensions is equalised prior to calculating said parameter 

30. A method of characterising ventricular operation of a patient's heart, comprising sensing a plurality of electrical 
signals associated with the propagation of a repoiarisatioh wave originating in the patient's heart from different 
spatial positions, with respect to the patient's heart., processing data obtained from the plurality of electrical signals 
using a decompositon technique to yield a plurality of vectors describing the repolarisation wave which are defined 
with respect to a set of axes in an optimised orthogonal domain space, one of the axes being aligned with respect 
to a direction of maximum energy for the repolarisation wave, wherein ventneuiar operation is characterised by 
calculating the proportion of repoiarisatioh energy not contained within a three dimensional space represented by 
the three most significant energy directions. 

31. A method of characterising ventncular operation of a patient's heart comprising sensing a plurality of electrical 
signals associated with the propagation of a repolarisation wave originating in the patient's heart from different 
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spatial positions with respect to the patient's heart, processing the plurality of electrical signals to yield a vector 
which is representative of the wavefront of the repolarisation wave with respect to a first set of axes,. transforming 
the vector to a second set of axes defining an optimised orthogonal domain having a first axis aligned with a 
direction of maximum energy, the domain comprising three dimensions representing the dipolar components of 
5 the repolarisation wavef ront vector and at. least one further dimension, representing the nqn<jipolar components 

of the repolarisation wavef rdnt vector r and determining the energy of the non-dipolar components. 

32. A method as claimed in claim 31 wherein said optimised orthogonal domain has eight dimensions and the trans- 
• . formed repolarisation wavefront vector S. has eight components, s^to sg cprresponding one to each dimension, 

io ... :..wherein ihevector component's are ranked in order of most significance with respect to energy and the non<iipqlar 
components are represented by the fourth to eighth components s 4 to s 8 . 

33. A method as cjaimed in claims 31 or 32, wherein the energy of the non-dipolar components is determined for a 
portion of the repolarisation wave corresponding to a particular region of the heart muscle 

75 

34. A method of determining whether a patient is suffering from heart failure, and preferably, a method of whether the 
patient is suffering from one of hypertropic cardiomyopathy, idiopathic dilated cardiomyopathy and acute myocar- 

. ; -dial infarction, comprising measuring the energy .of the non^ipolar components of a repolarisation wave produced 
by a patient's heart, wherein preferably this is accomplished by employing an implantable medical device which 
20 is capable of sensing a repolarisation wave and prpcessing.same. 

35. A method of determining depolarisation start and end points for measuring characteristics of a signal representing 
changes in energy during depolarisation. of a patient's heart, comprising finding a first peak in the energy of the 
signal corresponding to depolarisation ofthe patient's heart, determining a point in time, t RP , corresponding to the 

25 peak energy and determining the maximurrv energy E Rma i of the signal at that point, .determining a point in time 

t' RS before t RP and a point in time t' RE after tRp where the energy ofthe signal drops to a predetermined percentage 
of the maximum energy, determining the depolarisation start point by subtracting a first predetermined time interval 
from time f RS and determining the depolarisation end point, by adding a second predetermined time interval to 
timet' RE 

30 

36. A method as claimed in claim 35, wherein said predetermined percentage of the maximum energy E Rmax is in 
range of 50 to 90% of. E Rmax , preferably, in the range of 60 to 80% of E Rmax , and most preferably 70% of E Rmax . 

37. A method as claimed in claim 36, wherein said first and second predetermined.time intervals. are in the range of 
35 38 to 58 msec, preferably in the range of 43 lo 53 msec, and most preferably 48 msec. 

38. A method of determining repolarisation start and end pomts for a signal representing, changes in energy during 
depolarisation and repolarisation of a patient's heart, comprising: . 

40. finding a first peak in the energy of the signal corresponding to depolarisation of the patient's heart, determining 

the maximum energy E Rmax of the signal at the peak and determining a point in time t-' RE where the energy of 
signal has dropped, to a predetermined percentage of E Rrnax , finding the next peak in the signal energy cor- 
responding to repolarisation and determining the point in time tjp where that peak occurs, determining the 
repolarisation start point as a predetermined fraction of the time interval between t' RE and t TP , 

45 determining the repolarisation end point by determining, a vector s 2D (tj) which describes the repolarisation 

wavefront as projected oh to a plane spanned by two orthogonal vectors u T and u 2 :which represent the max- 
imum energy and next most energy of the repolarisation wave in two orthogonal directions for t, > t TS , the 
vector having a tip which defines a path on said plane, dividing the area defined by the path of the tip of s 2D 
(tj) in the plane of u 1 and u 2 into a plurality of equal rectangular cells, assigning a measure Dj lo each cell 

so dependent on the time spent by the tip of s^oft) in the 1 th cell, discarding cells having the measures E>, = 0 and 

ordering ail other cells in respect of D,, determining a threshold value D m of D r which is greater than the mean 
value of Dj, and determining the end point of repolarisation t TE as a point at which Dj > D, h 

39. A method as .claimed, in claim 38, wherein said predetermined fraction of the time interval between f RE and % is 
55 in the range 34 to Vi preferably said predetermined fraction of the time interval between f RE and t XP is 

40. A method as claimed in any of claims 1. to 9, or 20 to. 33, wherein the plurality of electrical signals are measured 
from electrodes positioned in different spatial positions corresponding to the standard ECG . leads of I, II, V1, V2, 
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V3. V4, V5 and V6. 

41 . A method as claimed in any of claims 1 to 9, 20 to 33 or 40, wherein data from the electrical signals is decomposed 
into an optimum domain space using singular value decomposition 

42. A method of detecting whether a patient is heaithy or sick, wherein more than one parameter characterising ven- 
tricular operation. as claimed in any ol claims 1 to 9, 21 to 33, 40 or 41 is calculated. 

43. ;A method of detecting abnormalities of ventricular repblarisation in a patient using a. method of characterising . 
ventricular operation as claimed in any of claims 1. to 9, 21 to 33; 40 or 41 ,.or.a method, of monitoring changes in 
the autonomic system of the patient as claimed in claim 20, -wherein preferably said method is a method of detecting 
the onset of ischemia. 

44. A method of determining whether a patient is suffering from heart failure, and preferably a method of determining 
whether the patient is suffering from one of hypertropic cardiomyopathy, idiopathic dilated cardiomyopathy and 
acute myocardial infarction, comprising measuring a parameter characterising ventricular operation as claimed in 
any of claims.'! to 9, 21 to 33,, 40 or 41 , wherein preferably the method is accomplished by employing an implantable 
medical device which is capable of sensing depolarisation and repolarisation waves, and processing same 

45. A method of categorising subjects using a method of characterising ventricular operation as claimed in any of 
claims 1 to 9, 21 to 33, 40 or 41, or using a method of monitoring changes in the autonomic system of a patient 
as claimed in, claim 20. 

46. A method of monitoring the autonomic system of a subject, wherein a method of characterising ventricular operation 
as claimed in any of claims 1 to 9, 21 to 33, 40 or 41 , is used to measure changes in the autonomic system of the 
subject, wherein preferab!y : the method is . accomplished by employing an implantable medical device which is 
capable of sensing depolarisation and repolarisation waves, and processing same. 

47. A method as claimed in claim 46, wherein changes in the autonomic system of the patient are measured to. 

a) monitor the progress of a disease in a subject; 

b) monitor the influence :of drugs oh the autonomic system of a subject, 

c) control the rate of pacing for a pacemaker; 

d) control a drug dispensing pump arid./ or 

e) test the response of the autonomic system of a. subject. 

48. A method as claimed in claim 47, wherein changes in the autonomic system of the patient are measured while the 
patient executes a pre-determined procedure to. effect a change in autonomic system, wherein preferably the 
autonomic system is measured while Valsalva manoeuvre is performed by the patient and/ or while a set of postural 
changes are executed by the patient, 

49. A method of detecting whether a patient is healthy or sick comprising sensing a plurality of electrical signals from 
different spatial positions in, on or near the patient's heart during depolarisation and repolarisation of the patient's 
heart and determining a parameter.. 

a) Which measures vector deviation of the depolarisation and. repolarisation waves in terms of Total Cosine 
RJOJT (TCRT), 

b) which measures the spatial variation of a repolarisation.wavef rpnt in terms of T~wave Morphology Dispersion 
(TMD), and / or 

c) which measures local repolarisation defects as a function of the energy of the nonrdipolar energy compo- 
nents of a repolarisation wave. 

50. A system for measuring characteristics of ventricular operation of a patient's heart in accordance with any of the 
methods of claims 1 to 9, 21 to 33, 40 or 41, or a system for monitoring changes in the autonomic system of a 
patient as claimed in claim .20, comprising a plurality of electrodes for sensing a plurality of electrical signals in 
different spatial positions in, on or near to a patient's heart to monitor the propagation of depolarisation and repo- 
larisation waves originating in the patient's heart, a microprocessor for processing the plurality of electrical signals 
to determine a characteristic of ventricular operation as claimed in any of claims 1 to 9, 21 to 33, 40 or 41, or 
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condition of the patient's autonomic system as claimed in claim 20, and an indicator to convey the value of the . 
measured characteristic, or condition of the patient's autonomic system, to an observer, wherein preferably the 
indicator is a visual display and/or an audible alarm 

51. A system as claimed in claim 5Q, wherein the system senses the plurality of electrical signals.via. 

a) at least one electrode attached to a medical electrical lead of an implantable medical device selected from 
the group of pacemaker, cardioverter, defibrillator, pacemaker^rdibverter-defibrillator and. heart monitor 
which is implanted within the subject's body; 

b) .at least 1 one electrode provided on a housing of an implantable medical device, 

c) at leasfone "sTibcLrtan^ 

d) at least one external electrode applied to the patient's body 

52. A system as claimed in any of claims 50 or 51 , wherein the system has eight electrodes which correspond to the 
standard electrode positions of I, II, VI, V2, V3, V4, V5 and V6 

53. An implantable medical device which measures characteristics of ventricular operation of a patient's heart in ac- 
cordance with any of the methods of claims 1 to 9, 21 to 33, 40 or 41 , or monitors changes in the patient's autonomic 
system as claimed in claim 20, comprising a plurality of electrodes for sensing a plurality of electrical, signals in 
different spatial positions in, on or near to a patient's heart to monitor the propagation of depolarisation and repo- 
larisation waves, originating in the patient's, heart, and a microprocessor for processing the plurality of electrical 
signals to measure a characteristic of ventricular operation as claimed in any of claims 1 to 9, 21 to 33, 40 or 41., 
or condition of the patient's autonomic system as claimed in claim 20, and to generate an output signal corre- 
sponding to; said measured, characteristic 

54. Apparatus as claimed in claim 53, wherein the apparatus is at least one of. 

a) a pacemaker or pacemaker-cardioverter-defibrillator having a pacing rate which is controlled in response 
to. the calculated characteristic of ventricular operation, . 

b) a cardioverter or.a defibrillator, the operation of which is controlled in response to the .calculated character- 
istic of ventricular operation, 

c) a monitor for detecting abnormalities in ventricular repolansation, the device preferably being able to trigger 
a warning signal when the calculated characteristic of ventricular operation is within certain limits, 

d) a monitor for measuring changes in the autonomic system of a subject, the device preferably being able to 
trigger a warning signal when the calculated characteristic of ventricular operation is within certain limits, and 
/or 

e) a device for controlling a drug dispensing pump 

55. A computer program, compnsing software code portions for performing the method as claimed in any of claims 1 
to 9, 20 to 33, 40 or 41 , a computer program product which is directly loadable into the internal memory of a digital 
computer, comprising software code portions for performing the method as claimed in any of claims 1 to : 9, 20 to 
33, 40 or 41 , a microprocessor which is programmed with software code portions for performing the method as 
claimed in any of claims 1 to 9, 20 to 33, 40 or 41 , or a digital computer which is programmed with software code 
portions for performing the method as claimed in any of claims 1 to 9, 20 to 33, 40 or 41 
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